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Abstract 


This position paper on Inadvertent Weather Modification is based 
on data characterizing the physical, chemical and dispersion state of 
the stabilized ground cloud vithin the first three hours after launch 
supplied by WASA-Lyndon B. Johnson Space Center. From this government- 
supplied information "we have estimated the extent of the S.G.C. and the 
particle content for times three hours, one day, three days and seven 
days after launch. We then discussed in detail uhe involvement of the 
S.G.C. in warm and cold cloud formation processes and the interaction 
of the S.G.C. aerosol with solar radiation for the times after launch 
mentioned above. Based on the climatology of the Florida Peninsula, 
we assessed the risk for weather modification. Certain weather situa- 
tions warrant launch rescheduling because of the risk of 

- possible impact on .hurricanes 

- hail formation and lightning activity 

- strong wind developments 

- intensification of high rainfall rates 

•The cumulative effects of 40 launches per year on weather modification 
were found to be insignificant. 
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Preface 


We have investigated the possible impact of the stabilized space shuttle 
exhaust cloud ("Stabilized Ground Cloud, S.G.C.") on the "weather" of 
the Florida Peninsula for a time period of three hours after launch up 
to seven days after launch. This position paper is based on information 
supplied by WASA-Lyndon B. Johnson Space Center (all data on the S.G.C.) 
and on information extracted from pertinent literature. An assessment 
team was formed consisting of the following members who have complemen- 
tary research experience in vital areas of inadvertent weather modification 


Dr. Volker A. Mohnen, Director 
Atmospheric Sciences Hesearch Center 
The University at Albany 


Dr. Vincent J. Schaefer 
Leading Professor, Atmospheric 
Sciences Research Center, The 
University at Albany 

Mr. Eugene BoUay, Former Chief, 
Office of Weather Modification, 
HOAA 

Dr. G. Garland Lala, Research 
Associate, Atmospheric Sciences 
Research Center, The University 
at Albany 

Dr. Patrick Souires, Head 
Atmospheric Physics, Desert 
Research Institute, Reno, Nevada 

Dr. James E. Jiusto, Head 
Atmospheric Physics, Atmospheric 
Sciences Research Center, The 
University at Albany 

Dr. Lance Bosart, Associate 
Professor, Dept, of Atmospheric 
Science, The University at Albany 


Scientific Project Director and 
Chairman of the Assessment Team. 
Principal discussant of Chapter VI 
entitled "Solar Attenuation Model 
for the Stabilized Ground Cloud." 

Principal discussant of Chapter I 
entitled "Historical Notes." 


Principal discussant of Chapter II 
entitled "Summary Review from. Overt 
Weather Modification Efforts . " 

Principal discussant of Chapter III 
entitled "Assimptions and Numerical 
Values." 


Principal discussant of Chapter IV 
entitled "Warm Clouds . " 


Principal discussant of Chapter V 
entitled "Cloud Physics Processes - 
Cold Clouds." 


Principal discussant of Chapter VII 
entitled "Florida Synoptic 
Climatology. " 



Dr. Earl G. Droessler, Reviewer 
ITorth Carolina State University 
Research Administration 
Box 5356 

Raleigh, North Carolina 


The assessment team iset twice for three days each at the Institute on 
Man & Science (April 1976) and at the Atmospheric Sciences Research 
Center's Field Station at Whiteface Mountain, New York (July 1976). 

In the interim period, the memhers of the team have "been in regular 
contact through individual and conference calls. The problem was 
approached as outlined in the block diagram. 
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The following risk situations for inadvertent weather modification 
due to the space shuttle exhaust were concluded: 

1. Exhaust cloud encountering active convective precipitation cells 
with consequent vertical transport to the upper troposphere and potential 
for acid rain 

(a) sea hreeze convergence during the warm season with attendent 
afternoon thunderstorms. Effects include possible localized 
hail and brief wind gusts in excess of 20 ms~^ . Affected area 
is less than 100 km^ with a time scale of less than T + 1 day. 

(b) frontal and prefrontal activity including squall lines -witli 
attendent thimderstoimis. Effects include possible localized 
hail, wind gusts in excess of 30 ms~^ and tornadoes. Affected 
area is 100-500 km^ with a time scale of less than T + 2 days. 

(c) general air mass thunderstorms not associated with (a) and (b) 
above but responding to different s\mmer synoptic flow patterns. 
Effects include possible localized hail and brief wind gusts 

in excess of 20 ms~^. Affected area is less than 100 km^ with 
a time scale of less than T + 1 day. 

(d) tropical storms in the vicinity of the Florida peninsula 
within 2k hours of launch time. Potential effect of shuttle 
exhaust cloud caught up in the circulation of a tropical storm 
is unknown in terms of inadvertent weather modification.' A sub- 
sequent change of direction of such a storm might be interpreted 

aJi "act of God" by some people with possible social and 
legal problem from communities in the landfall region. 

2. In the months Uovember-April , when advective and radiative fogs 
maximize, very significant worsening of visibility conditions in foggy 
situations could occur within the area affected by the dissipating S.G.C. 
up to T + 1 day (area affected up to 10^'' km^) and particularly under wind 
flow conditions from the S-E quadrant. 

3. Minor risk associated with easterly flow in lower troposphere 
(unless tropical disturbances are present) particularly in those situations 
where atmosphere is table- under those conditions, clouds do not reach 

the level where ice phase processes are operative. However, overseeding 
of warm clouds with CCK could result in a very significant reduction of 
precipitation over the entire area affected by the dispersing cloud. 

Effect diminishes after T +_ 1 day, (Criteria: shallow warm cloud 

system and no ice phase.) 

U, Stagnating anticyclonic conditions with reduced dispersion of 
S.G.C. Little cloudiness is normally associated with conditions of this 
type. The impact is therefore restricted only in the area of visibility 
deterioraticnand solar energy reduction. This constitutes therefore a 
nuisance and conceivably might violate EPA standards. On rare occasions 
air mass thunderstorms may develop, particularly along the sea breeze 
convergence zone, under stagnant anticyclonic conditions during the 
warm season. The risk would then be equivalent to l(c) above. 


iii 



5. Possible modification of a major hurricane located east of 
Florida peninsula at time of launch. Air from launch site -would participate 
in the storm circulation and might indeed cause some modification effects 
producing unkno-wn results. Any subsequent veering of such a storm -would 
undoubtedly cause serious social and legal problems. 

6. Cumulative effects: for the projected 40 launches per year 

assuming several days spacing bet-ween laimches is considered negligible. 

7. Minimal risk and impact: strong -westerly -winds system extending 

through the lower troposphere. 

Certain -weather situations warrant launch rescheduling because of the 
risk of 

- possible impact on hurricanes 

“ hail -formation and lightning activity 
~ strong -wind developments 

- intensification of high rainfall rates 

The cumulative effects of launches per year on -weather modification 
■were found to be insignificant . 
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HISTORICAL NOTES 

The impact of modern technology on the quality and properties of 
the global atmosphere has been the continuing interest of investigators 
for more than 35 years. 

Some of the first problems studied related to the role that small 
airborne particulates played in the formation of snow and ice crystals 
in the atmosphere of eastern New York (Schaefer, 1942). This expanded 
within a few years to the subject of the filtration of poisonous smokes 
(Langmuir and Blodgett, 19^1 ) and then the utilization of large 
quantities of smoke to protect army troops in Europe and navy ships in 
the South Pacific (Langmuir, 1948). 

There followed investigation of the role of airborne snow on the 
communication effectiveness of airborne radios on B-IJ and other air- 
craft (Schaefer, 1947), then of the safety of aircraft flying in super- 
cooled clouds (Langmuir, et al. , 1946), along with the study of snow 
crystals and methods for replicating them in the free atmosphere 
(Schaefer, 1941) which led to the discovery of a method for the seeding 
of supercooled clouds in the free atmosphere (Schaefer, 1946, 1950). 

During the period of 1944-45 some studies were carried out on the 
effectiveness of AI 2 O 2 as an ice-forming nucleus, and the role played 
by organic and inorganic vapors on the growth habits of ice crystals 
(Schaefer, 1953). 

The possible role that particles and gases play in local, regional 
and global weather patterns is a research concern of an increasing 
amber of scientists. Some of the work was summarized in a review paper 
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several years ago (Schaefer, I 969 ). More recently it vas the theme of 
a major conference on inadvertent weather modification (Blanchard, 19T^). 

Studies were conducted on the effectiveness of several. different 
types of particles as heterogeneous ice nuclei (Schaefer, 19^9; 

Vonnegut, 19^7) which followed shortly after Schaefer's discovery 
that dry ice and other substances colder than -40°C(F) led to the forma- 
tion of tremendous numbers of ice embryos formed throxxgh homogeneous 
nucleation. Shortly after these discoveries were announced and sub- 
stantiated by field work. (Schaefer, 19^7) in supercooled clouds, Langmuir 
proposed that warm clouds could also be modified through the sudden 
injection of large water drops or hygroscopic materials that would trigger 
or enhance coalescence and thus initiate the precipitation process 
(Langmuir, 19^8). At a considerably earlier time, Houghton and Radford 
( 1938 ) had conducted field studies of the stability of ground fogs. 
Following Schaefer and Vonnegut 's work with dry ice and silver iodide, 
Fournier d'Albe, et al. (1955) in India carried out some seeding of warm 
clouds using sodium chloride. Langmuir had suggested the use of ammonium 
sulfate for similar purposes , but the field tests we carried out were 
not very encouraging. 

Since that time, others have attempted the modification of warm 
clouds by salt particle seeding (Schleusener and Kocielski, 19715 
Biswas and Dennis, 19715 Simpson and Dennis, 197^)- It shoinld be noted 
that the utilization of hygroscopic materials to affect the precipitation 
process in warm clouds either overtly or inadvertently, involves an 
entirely different mechanism than is involved when a change of phase 


is involved. 
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The Tise of salt-like or watery substances involves particles that 
must be larger than 25 pi to be effective. Since dry ice or silver 
iodide seeding in supercooled clouds can be effective with particles less 
than 0,01 pa, this means that a difference in' weight or mass of about 
10 billionfold is involved in producing effective particles in the two 
processes if one assumes that there is a one to one relationship between 
nucleus and precipitation particle. Langmuir has suggested that the 
injection of large water droplets into an actively developing warm cloud 
would lead to a chain reaction (Langmuir, 19^8). Thus, a relatively 
small amount of seeding material could dramatically accer.tuate the 
growth mechanism. More recent attempts to model such growth mechanisms 
suggest that this may be possible (Berry, I96T). 

In any event, it is important for this report on the space shuttle 
program to consider both effects, i,e., whether either or both warm and 
cold cloud modification is possible from the residue of the ground cloud 
produced during launching. 

The concentration of airborne particulates affecting the air quality 
of the global atmosphere has been a recent study of one of us (Schaefer, 
1972; 1973; undated). Earlier it was pointed out evidence exists that 
the effluent of sugar centrales in Puerto Rico produced (Schaefer, 

19^9) much larger convective clouds than formed over the sea upwind of 
the island under tradewind- cloud development. It was also suggested 
(Schaefer, 1958) that the massive convective clouds which formed in 
Africa prior to the start of the rain season showed strong evidence of 
being overs eeded with cloud condensation nuclei produced by the burning 
of the bush on a massive scale. Such clouds, though forming to vertical 



heights of 35-40,000 feet, produced little if any precipitation, thou^ 
the tops of these clouds formed ice crystal streamers extending for 
hundreds of miles. 

The possibility of the overseeding of clouds by the launch cloud 
products, ■which might enhance their stability and repress the preci- 
pitation-forming mechanisms in the local and downwind atmosphere, ■will 
also be considered. 

Figure I-l shows the basic relationships which exist between air- 
borne particulates on the global scale as related to natural and 
anthropogenic sources. Figure 1-2 then shows typical patterns found in 
maritime air, and Fig. 1-3 in continental air. Since the air flow 

pattern at the lower levels of the atmosphere over the Florida peninsula 

\ 

are likely to consist of maritime air, i'ts modification toward continental 
air properties would likely be of interest in assessing possible 


inadvertent effects 
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TITLES 
Figure I-l: 
Figure 1-2: 
Figure 1-3: 
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The Basic Relationships of Particle Sizes which Exist in 
the Free Glohal Atmosphere 

Typical Particulate Sizes and Concentration Found in 
Maritime Air 

Typical Particulate Sizes and Concentration Found in 
Continental Air 
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Chapter II 


SUMMARY REVIEW 

FROM OVERT WEATHER MODIFICATION EFFORTS 


A, BACKGROUND 

Soon after the discovery "by Langmuir, Schaefer (19^6) and Vonnegut 
( 1947 ) that there "was a scientific basis for man to modify nature’s way 
of producing rainfall, economic pressirres and entrepreneurs were ready 
to put this knowledge to use. The time was 1946. ^“Jbere are we now, 
thirty yesirs later? 

Weather modification is an all encompassing term, still way beyond 
man’s fondest dreams. What we can do is cloud seeding. It is a 
reality and there is scientific theory to back it up. 

The underlying physical basis for cloud seeding is the altering of 
the size spectrum or phase of the cloud condensate through the manipula- 
tion of the popilLation of the condensation, freezing or sublimation nuclei. 
There is experimental evidence that micro physical changes and dynamic 
changes are brought about when nuclei are added or when naturally 
occurring nuclei are activated. 

The internal characteristics of clouds determines the behavior of 
clouds, whether they grow, whether precipitation elements can develop 
from the cloud droplets, and the lifetime of the clouds. The nature of 
the ■ aerosols and the updraft speed determines the number of cloud droplets 
formed with the aerosols serving as condensation nuclei. Different air 
masses, maritime and continental, have varying aerosol populations. 
Typically, maritime air laasses have aerosols ranging from a few tens to 
hundreds per cm^ , whereas continental air masses may vary from several 
hundred to several thousand per cm^. The mass of average cloud droplets 
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is about 10”^ gram and the mass of an average precipitation element is 
about 10“^ gram. 

There are two processes by means of which we can explain the trans- 
foimation and growth of cloud particles into precipitation drops. 

One process leading to the formation of precipitation elements 
depends on the availability of ice particles and the presence of super- 
cooled cloud particles. Because ice has a lower vapor pressure than 
water there is a transfer of water vapor to the ice particles leading 
to the growth of precipitation size elements . This is known as the 
Bergeron-Findeisen process. It is believed to be the principal mechanism 
in causing precipitation in the middle and high latitudes. 

Another process explains the formation of precipitation elements 
from clouds which are warmer than freezing. The process depends on the 
collision and coalescence of cloud droplets having different size and 
mass falling at different speeds through the air. In maritime air> 
particularly where larger cloud droplets are expected because of the 
presence of giant condensation nuclei » the larger cloud droplets can 
grow at the expense of the small droplets very quickly to reach raindrop 
size and fall out as precipitation. In air mass clouds of continental 
origin there are many more condensation nuclei and there is very active 
competition for available moisture amongst a large population of cloud 
dropleus and frequently no precipitation sized particles form and, 
therefore, iw rain falls from many such clouds. 

It is believed that both processes, the Bergeron and the collision 
and coalescence process are active- and in competition simultaneously. 
There are no adequate statistics that establish the percentage of the 
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total precipitation produced by either process or -whether there is a 
critical time factor -when one process may be the dominant one in the 
life cycle of a cloud. 

An additional feature of the Bergeron-Findeisen process leads, to 
an increase in temperature resulting in increased cloud buoyancy vhen 
supercooled cloud drops are converted to ice crystals and the latent 
heat of fusion is released. 

The impact man can make in -weather modification is to affect either 
the condensation nuclei population and the ice forming nuclei population 
in clouds (i.e., aerosols). Adding nuclei of one kind or another is 
called cloud seeding. 

There is ample evidence that precipitation -^fill decrease when there is 
an overabundance of nuclei, either cloud forming nuclei (CCIT) or ice forming 
nuclei (IN) (Project White Top, Climax I, etc.). Almost all cloud seeding 
experiments in the past were carried out in the anticipation that htere was 
a shortage in natural nuclei which is not always the case. In California 
and Colorado where cloud seeding has been undeirray during many winter seasons , 
clouds ha-ving optimum seedability have been discovered as a function of 
cloud top temperature (Grant and Elliott, 197^) • VJhenever clouds are seeded 
in storms ha-ving sufficient depth to reach a cloud top temperature in the 
range of -5°C to -25°C, increases in precipitation occur. When colder 
cloud top ten5>eratures are observed and the clouds are seeded, decreases 
in precipitation will fodlLow, suggesting, that there are .so many natural 
ice nuclei that the addition of artificial nuclei will only increase the 
coE5>etition for the available water vapor and fewer ice cs^rstals can grow 
to the size of precipitation elements. Figure II-l illustrates 
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schematically the cloud top ten^erature window. 

Figure II-l CLOUD TOP TEMPERATUEE 
SEEDING NXiJDOW 

PRECIPITATION 
RATE 


The most common nuclei used for cloud seeding are Salt and Urea for 
warm cloud seeding (Collision and Coalescence Process), and Silver 
Iodide for cold cloud seeding (Bergeron-Findeisen process). 

Precipitation Modification using the warm cloud seeding process 
(Collision and Coalescence) has received relatively little attention in 
this country and in the rest of the world. There are no long time 
experiments that have been accepted by the scientific community as 
yardsticks of performance or achievement. Research experiments of 
short duration and operational programs have been carried out , some 
with encouraging results but statistically inconclusive. 

One reported experiment (Biswas, et al., I967) using salt as a 
warm cloud seeding agent was conducted for a total of I8 ejip>eriiaental 
seasons in three climatologically similar regions in northwest India during 
the years 1957 - 1965 . In this statistically randomized experiment, a 21 ^ in- 
crease in the total seasonal precipitation was reported. 
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The technical feasibility of modifying vrarm fog has been demon- 
strated in certain meteorological situations. Results have shovn that 
the application of heat to disperse -warm fog at airports appears to be 
a more direct and dependable technique without undesireable side effects, 
such as corrosion from salt. 
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B. MHEES DO WE STMD IH COLD CLOUD SEEDING 

1. Cold Fog 

Cold fog dispersal is operational at selected airports and air 
bases using botb surface and aircraft type silver iodide dispensers. 

2. Snowpack Augmentation 

SnoTTpack can be increased and in some cases redistributed in 
mountainovis areas of the vest by seeding vinter orographic clouds. One 
of the longest cloud seeding programs (Elliott, 19T5) carried out in 
this country, and still ongoing, shovs an 8 * 5 % increase in annual stream 
flov for over a 25 year period. The seeding has been carried out largely 
vith ground generators dispersing silver iodide particles in the size 
range of 0.01 to 0.1 microns diameter. During a typical operational day,\ 
covering 15 hours, about 630 gr. of silver iodide are released from a 
total of seven generators scattered throughout the mountains. 

3. Mesoscale Modification (Winter Convective Band Modification) 

A program (Brovn and Elliott, 1972) indicating that mesoscale 

changes can be produced by seeding convective bands imbedded in vinter 
storm clouds vas carried out in Santa Barbara county during 1967-1971* 
Statistical evidence suggested that precipitation increases over 100^ 
vere occurring over large areas 100-150 Jniles dovnvind from the seeding 
source. The seeding nuclei vere produced vith LW-83 pyrotechnic devices 
(Navy) producing 400 grams of silver iodide in a three minute period. 

These devices vere burned every 15 minutes as long as the convective band 
vas over the seeding area. This program can be considered a forerunner to 
programs related to the modification of vinter storm systems on the Pacific 
Coast. It also presents good 'evidence that sizeable extra-area effects 


are likely and possible. 
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4. Cumulus Modification 

One of ' the classic experiments (Simpson and Woodley, 1971) in 
weather modification deals with the precipitation augmentation from 
carefully selected individual cumulus clouds in south Florida. Seeding 
effects of 200 to 300^ have "been demonstrated. The seeding operations 
were carried out from aircraft and the silver iodide nuclei were generated 
hy pyrotechnic devices. The amount of silver iodide expended per day of 
experimentation typically ranged from 4 to 15 kilograms. The same 
program is now investigating the precipitation augmentation from 
multiple convective clouds and assessing the total area effect . The 
work is still in progress. Preliminary indications suggest much smaller 
increases in precipitation by seeding multiple clouds. 

The study also showed that in south Florida (Holle, 19T4), there 
are cumulus clouds favorable for seeding from mid-April through mid- 
September. Outside of this season seedable cumulus clouds are more 
rare and at least half of them are associated with cold frontal activity 
moving down the peninsula. 

5. Hail Suppression 

Commercial operations related to hail suppression activities report 
success but lack satisfactory evaluation systems. A major National Hail 
Research Experiment has been established by the National Science Foundation 
at the National. Center for Atmospheric Research. As of this date no 
results can be reported. Modification of hail forming clouds is reported 
to be operational in the U.S.S.R. and to some extent in South Africa. 
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6. Hurricane Modification 

THe modification of hurricanes (Gentry, 1970 ) hy expanding the eye 
of the storm radially outwards by cloud seeding is one of the major 
national goals . Experiments in hurricane Debbie in 19^9 provided en- 
couraging evidence that maximum winds can he reduced hy a pattern of 
cloud seeding using silver iodide. Plans are now being made and the 
necessary eq^uipment and instrumentation procured to carry out an 
important demonstration program. 

7. Lightning Suppression 

Lightning suppression experiments (Kasemir, 1973) using silver 
iodide are inconclusive. Another method of lightning suppression still 
in the experimental stage uses "chaff needles" to bleed off the elec- 
trical fields. This appears a promising method if operational means 
for delivering the chaff to the electric field forming centers in 
cumulo-nimbus clouds can be developed. 

The modification of tornadoes or other severe storms are research 
interests still in the design stage. Some excellent modelling efforts 
are underway to deteimine where and when severe storm modification 


might 'be feasible. 
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Chapter III 


ASSUMP.TIOJIS Am) irai^SRICAL VALUES 


A. IKIEODUCTIOIT 

In any attempt to assess the weather modification effect of an 
aerosol cloud, one of the most important aspects is the particle size 
distribution and number concentration. The distribution and concentra- 
tion are important in determining whether entrainment of the aerosol 
into a cloud will inhibit or promote precipitation and will also 
determine the dominant precipitation mechanism. Equally important is the 
background aerosol character and its concentration relative to the’ 
aerosol being introduced. Knowing the distribution and concentration of 
the two aerosols allows one to estimate whether there will be any 
weather modification impact and the magnitude of the effect. 
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B. BACKGEOOTD AEROSOLS 

From a ■warm-cloud weatlier modification perspective j specification 

of the aerosol in terms of a supersaturation spectrum is more meaningfiol. 

A supersat\aration spectrum gives the total number of particles (cloud 

condensation nuclei) activated at a given supersaturatlon. Typical 

supersaturations used are in the range of 0.2^ to 2$ -which covers the 

range of supersaturation occurring naturally in clouds and fog. 

Supersaturation spectra follow a power law of the form 

fC 

N = CS 

where S is the supersaturation, C is the concentration at super- 
saturation and K is the slope of the spectrum. Measurements of the 
supersatiirate spectra over the Florida peninsula have been carried out 
by Fitzgerald (1972) for aerosols of both maritime and continental 
origins. The res'ults of these meas'urements are s-ummarized below 
Maritime aerosol N = 2h20 

Continental aerosol R = 5913 S°‘^^ 

where R is the concentration of particles (cm“^) activated at a super- 
saturation S, the latter being expressed in absolute units; for example, ’ 
at S = 10”^ (often expressed as ”1^") R = 291 cm“^ in a maritime 
aerosol, or 515 cm“^ in a continental aerosol. 

Backgro-und levels of ice nuclei for evaluating the potential 
weather modification impact on cold clouds are given in Chapter V 
(Figure V-2). 
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C, VOLUME OE THE SGC ASH) AEROSOL MASS CONCEi3TRATIOH 

1. Cloud Volume as a Function of" Time 

The specification of the volume of the stabilized ground cloud is 
based on a linear extrapolation of RASA measxirements of the cloud volume 
shortly after launch (NASA Staff -Langley Research Center). Volume 
estimates derived in this manner are probably conservatively smal 1 but 
offer a reasonable approximation considering the lack of data on the cloud 
volume at long times after the launch. Note that extrapolation of the growth 
curve plotted along with the data leads to a cloud volume of 1200 km^ at 
3 hours as opposed to the volume of 600 km^ used in these calculations. 

Table III-I below gives the estimated volume of the stabilized ground cloud 
at various times after launch. 


TABLE III-I 

Estimated Cloud Volume vs. Time 


Time 

T + 3 hrs. 
T + 1 day 
T + 3 days 
T + T days 


Volume of the Cloud 
6 X 10^ km^ 

5 X 10^ km^ 

1.5 X 10^ km^ 

3.5 X 10*^- km3 


2. Aerosol Mass Concentration in the Stabilized Groiond Cloud 
Specification of mass concentrations in the stabilized ground cloud 
is a difficTilt problem because of a lack of data on cloud volume and mass 
concentration at long times after the launch. The approach followed here 
is to start with the total emission of AlgO^ in the ground cloud and to assume 
that this mass is maintained in the cloud volume as specified in Table III-I. 
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The determnation of the total mass of Al^O^ in the stabilized 
ground cloud ■yr&s obtained from integration of the ceilculated emission, 
of the space shuttle vehicle as given in NASA-JPL Tech. Memo. 33-712. 
Based on the information that the material emitted during the first 20 
to Zh seconds (2 ioa altitude increment) makes up the stabilized ground 
cloud leads to an estimate of 10® grains of AJ. 2 O 3 in the cloud. Table 
III-II gives the mass concentration as a function of time based on the 
conservation of the total mass of material in the ground cloud and the 
cloud volume as given in Table III-I. (A 10 percent mass loss due to 
surface deposition is allowed during the first 3 hours . ) 

TABLE III-II 

AI 2 O 3 Mass Concentration vs. Time 

Time Mass Concentration of AI 2 O 2 

T + 3' hrs. 150 pg m“® 

T + 1 day I 8 pg m~® 

6 pg m“® 

2.6 pg m“® 


T + 3 days 
T + 7 days 



III-5 


D. AEROSOL SIZE DISTRIBUTION AND NUI4BER CONCENTRATION 

1. Size Distribution. 

Aerosol size distributions for the AI 2 O 2 aerosol were extracted from 
the data presented by Varsi (19T6). This data, presented in the form of 
histograms, was normalized and integrated to obtain cumulatiTre distri- 
butions of the fraction of particles larger than the indicated size (diameter) 
and is presented in Figures III-l through III-6. Power law functions were 
fit to these distributions with the resulting values for N, the fraction 
greater than diameter, in terms of D, diameter in microns, indicated on the figures. 

Inspection of these figures indicate that the data can be partitioned 
into two groups: one with a steep slope (Figoires III-l, III-3, III-5) 

associated with early measurements and the other with a shallow slope 
(Figures III-2, Ill-ii-, III-6) characteristic of later measurements or 
long time averages. This change of slope with time reflects the effects 
of coagulation at high concentration resulting in fewer small particles. 

For the purpose of this study, which is concerned with the cloud at long 
times after the launch, the distributions with the shallow slopes are 
more appropriate. 

Distributions obtained with the electrical mobility analyzer 
(Fig. III-3 through III-6) exhibit a shallow slope in the size region 
less than O.OT ym. This is typical of size distribution at small sizes 
where coagulation rapidly acts to remove the very small particles. 

Thus, to realistically describe the characteristics of the AI 2 O 2 
aerosol two distributions are needed; one for the small sizes 
(0.01 ym < D < 0.07 ym) and a second for the large sizes 
(0.07 ym < D < 50 ym) . The distribution at large sizes has been 
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Figure III-2 

Cumulative normalized data from wire impactor 








TION GREATER THAN DIAMETER 


III-IO 


-3.93 

N=L+.20E -5D 



DIAMETER CM1CROMS1 


Figure IH-5 

emulative nomalized data from electrical mobility analyzer (ETR- T + 7) 
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Figure III-6 

Cumulative normalized data from electrical motility analyzer (ETR- T + 13) 
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terminated at 50 Jim TDecau.se sedimentation of particles larger than this 
size •will rapidly remove them from the cloud. 

Inspection of the size distrihution at longer times (Figures III-2, 
III-4, III-6) indicate that an average distrilDution of the form 
~aw = D“^*^ dD (0.07 pm < D < 50 pm) 

is appropriate. At the smaller sizes , the data can "be described by a 
similar function of the form 

-dW = Uq dD (0.02 pm < D < 0.07 pm) 

where the slope is one half the value for the distribution at larger 
sizes. • 

2. Aerosol number concentration 

Having specified the form of the size distribution, the problem of 
determining the number concentration remains. The determination of the 
number concentration is best accomplished through the reguirement of 
mass conservation. By integration of the mass represented by the two 
distributions, assuming the density of an aerosol particle to be 
2.5 g cm~~ (Varsi, 1976), and the reguirement of continuity of the two 
cumulative distributions , the number concentration can be calculated 
from the mass concentration as given in Table III-H, Following this 
procedure , the number concentration was calculated as a function of time 
with the results given in Table 
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TAELS III-III 

Computed -Aerosol Concentrations (0.02 ym < D < 50 ym) 


Time 

T + 3 hrs. 
T + 1 day 
T + 3 days 
T + 7 days 


Aerosol JTuniber Concentration 

6.63 X 10^ m~^ 

7.96 X 10 ^ 

2.65 X loS 

l.li^ X 10^ m~3 ' 


Sor the purposes of this evaluation, the total mass of material in 
the cloud at 3 hours is assumed conserved at later times and the form 
of the distribution function is assumed constant with time. 
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S, SUMMARY AUD EECOMMSi®ATIONS 

Specification of the AI2O2 aerosol in the stabilized ground cloud 
has been accomplished in terms of cloud vol\ime, average mass concen- 

i 

tration, and the particle size distribution. These results have been 
derived from extrapolated cloud volume data, mass of AlgO^ 
ground cloud, and measured size distributions. Under the assur^tion 
that all mass present at three hours is present at later times and the 
size distribution remains constant in form, the resulting set of 
specifications are consistent with mass conservation. 

Future measurements of the aerosol in the stabilized groxind cloud 
should concentrate on the problem of. determining the mass balance 
and size distribution for times extending to at least several hours 
after launch. It would ad.so be of vital importance to meastire total 
particle concentration, and particularly the spectra of cloud conden- 
sation nuclei and ice nuclei during the same period. The availability 
of a complete and consistent set of data on uhe aerosol in the stabilized 
ground cloud is essential to improving estimates of the inadvertent 
weather modification impact of shuttle exhaust products. Measurements 
of size distributions aTone would be a quite inadequate substitute for 
specific measurements of those properties which are of direct cloud 


physical importance. 
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Chapter IV 


WARM CLOUDS 


A. IWTRODUCTORY REMARKS 

As described “in' Chapter' III, at 3 hours from launch, the 
SGC has a voliome of about 600 km^, having by then mixed with 
many times its original voliime of ambient air within the surface 
mixing layer. The enlarged and diluted cloud will then contain 
an aerosol consisting of the natural background aerosol, with 
the addition of particles derived from the SGC, the properties 
of which have been discussed in Chapter III. The 3-hour cloud 
will continue to be further diluted with additional volumes of 
ambient air during its subsequent drift downwind, and at some 
point in its history m.ay become involved in upcurrents which could 
result in the - formation of cloud and precipitation. Even while 
still in the surface layers, the SGC may become involved in con- 
densation processes if fog forms. 

The microphysical properties of the resulting water cloud 
may be modified in several ways, as a result of the presence of 
the artifically added aerosol. In some circumstances, such modi- 
fication could have significant effects on the formation of pre- 
cipitation and on the dynamical behaviour of the cloud. 

With the Florida area, the enlarged and diluted SGC v;ill usually 
become involved in convective (rather than stratiform) clouds. 

These form in updrafts caused by heating contrasts or by squall 
line convergence. The artificial aerosol can then influence cloud 
physics processes in the following ways; 
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(a) It may provide giant hygroscopic particles which act 
as the embryos for raindrop formation; 

(b) It may modify the basic micro structure of the cloud, 
tending to convert clouds which naturally would have had a 
colloidally unstable maritime- microstructure into colloidally 
stable continental type clouds, thus tending to delay rain forma- 
tion; 

(c) It may catalyze the formation of ice particles, so 
enhancing precipitation formation (discussed in Chapter V) ; 

(d) As a result of increased ice formation, the additional 
latent heat released may warm the cloud sufficiently to cause it 
to break through a stable layer.. If this occurs, the basic up- 
draft may increase, and the cloud may become both deeper and wider 
than previously, with the possibility of greatly increased pre- 
cipitation (Chapter V) . 
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B. THE S.G.C. PARTICLE SIZE DISTRIBUTION AND SURFACE PROPERTIES 
1. Size Distribution 

As discussed in Chapter III/ measurements made on the S.G.C. 
aerosol indicate that in the size range 0.035 ym<r<5 ym, the concen- 
tration {'K(rp)) of particles per cm^ for which is given by 

a relation of the form- 

In N(r ) = a - 2.5 In r 
Pi P 

It is assumed here that this same law holds for sizes up to 

X - 25 ym. In the region below r = 0.035 ym/ it is expected that ' 

coagulation in the early SGC will have reduced the slope of the 

distribution, which will be approximated dov?n to r = 0.01 ym by: 

In N(r_) = a^ - 0.75 In r 

p o P • - 

These_relations obviously must give the same value of N(r ) 

P 

when r = r^ = 0.035 ym. Thus, the constants a, and a , which 
are functions of time and decrease as the cloud is mixed and 
diluted, are related as follows: 

exp a^ = r^ exp a^^ 

= 3.531x10^ exp a^^ , 

all lengths being expressed in centimeters. 

Writing A(t) for exp a^, it then results that in the' 
range 0.035 ym<rp<25 ym, 

N(rp) = A(t) 

while in the range 0.01 ym<r <0.035 ym, 

P 

N(r^) = 3.531x10^ A(t) r 
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Taking the density of the aerosol particles (Chapter III) as 
-3 

2.5' g cm for all sizes, and writing for the mass concentra- 

< —3 ’ 

tion (g cm ) of the aerosol in the form of particles for which 
in the range 0.035 ijm<rp<25 ym the mass distribution is 
given by: 

M^(r^) - M^(r2) = 52,36 A(t) (r2°*^ - 

— 3 

In particular, the total mass of concentration (g cm ) in the 
whole range considered, from r = 0.035 ym to r = 25 ym, is 2,520 
A(t) g cm 

In writing similary M^(rp) for the range below 0,035 ym, 

^^o^^l^ “ %^^2^ ^ 1.233x10^® A(t) 

and over the range 0,01 ym<rp<0.035 ym, the total mass concentration 
is 6.143x10 ^ A(t]| . Thus the total mass concentration of the 

aerosol^_in the range r^ = 0.01 ym to 25 ym is 2.526 A(f) g cm"^. 

This expression for the total aerosol mass concentration is essen— ■ 
:ially independent of the lower limit of particle radius. 

At T t 3 hours, as discussed in Chapter III, the total mass 
of the aerosol particles in the SGC cloud is taken to be' 9x10^ g, 
the cloud volume at that time being 600 km^, so that the mean 
mass concentratxon is 1.5x10 g cm , Thus at T + 3 hours, 

A(t) = 5,938x10 . At later epochs, A(t) will be smaller, in 

inverse proportion to the total volume of the mixed and diluted 


SGC, 
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Thus, at T + 3 hours, in the size range 0.035 um<r <25 ym, 

N(r ) = 5.938x10“^^ r 
P P 

and over this range, the total number concentration is 2.591x10 
cm”^. 


In the range 0.01 ym<r <0.035 ym, 

P 

N(r ) = 2.097x10“^ r 
P P 


3 

and over this range, the total number concentration is 4.040x10 

-3 

cm . Thus the total number concentration at T + 3 hours in the 

3 3 

range 0.01 ym<r^<25 ym is 6.631x10 particles per cm . 

2 . Surface Properties 

It is assumed here that the surface properties of the AI 2 O 2 
particles have been irreversibly modified by reaction with gaseous 
HCl and H 2 O to form a soluble surface layer- The reaction is sup- 
posed to proceed to a uniform depth x cm in every particle, irre- 
spective of its size. It is further supposed that the reaction 
combines y grams of gaseous material with each gram of AI 2 O 2 / 


resulting in the formation of (1+y) grams of soluble hydrated sur- 
face chloride, which is assumed to have a constant effective mole- 
cular weight {in dilute solution) of M. Nominal values assumed 
— 8 "7 

are x = 10 and 10 cm, y = 1, M = 30. These values imply an 

areal surface density of the soluble product of 5xl0“^ to 5x10“”^ g cm” 

2 

(0.5 to 5 mg per m ) . 
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C. G-IMT AEROSOL PAETICLEB 

When a giant hygroscopic nucleus {r>lym) is exposed to the 
saturated air within a cloud, it forms a droplet which out- 
grows its neighbors which typically contain nuclei with radii of 
less than 0.1 ym; this is a result of the strong depression of 
the equilibrium vapor pressure over -the surface of the large drop- 
let. Such fast-growing droplets can become the embryos of rain- 
drops, since some of them can rather quickly reach sizes at which 
collision-coalescence growth begins. It is generally agreed that 
those giant soluble particles (such as sea salt particles) which 

have dry radii exceeding 10 ym are candidates to form precipitation 

-3 

embryos, and that if they occur in concentrations exceeding 10 

-3 .... 

cm , they may well cause precipitation formation. 

As discussed in paragraph B, the AI 2 O 2 particles derived 
from the SGC are only partly soluble.' Since the giant particle 
effect depends on the strong and persistent vapor pressure lower- 
ing by the solute (which is very large compared with the opposing 
Kelvin effect) , the relevant factor is the mass of soluble material 
in the particle. 

It is assumed here that an outer layer of depth x cm on each 
particle has combined with gases in the SGC to form a soluble 
material, the reaction involving y grams of gas for each gram of 

Al202- The resulting layer of soluble material then has a mass 

2 -3 

4ttp (l+y)r X, where, p“ is t^ken as 2.5 g cm (Chapter IV). 

XT tr 



IV-7 


-3 

The mass of a sea salt' particle of radius 10 cm and den- 
-3 -9 

sity 1.2 g cm is about 5x10 g. The modified Al^O^ particle 

which has a coating of soluble material of equal mass must have 

a radius of about (5x10 ^/47rp (l+y)x)^'^^. Taking P •= 2.5 g cm 

P P 

and giving y the nominal value of 1, and x a- maximum value of 

-7 . . -2 

10 cm, this is about 2.8x10- cm, which lies well above the assumed 

""3 

maximiam size (r = 2.5x10 cm) of particles in the SGC at T + 3 

IT 

hours. Even if the size distribution assumed for particle radii 

. -2 
exceeding r^ = 0.035 ym were extrapolated to 2.8x10 cm, the 

— 7 -3 

expected concentration at T t 3 hours would be only about 5x10 cm , 

which is very small compared with the concentration of such par- 

-3 -3 

tides necessary to initiate rain formation ("ulO cm ) . 

It may be concluded that at and beyond T + 3 hours, any giant, 
particles present would have no significant effect on rain forma- 


tion. 
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D, CLOUD FOEMATIOU 

The most critical property of a particle with regard to cloud 
formation is its critical supersaturation (S^) , that is, the maxi- 
mum supersaturation with which a large droplet can be in equilibrium. 
As discussed in paragraph b, it is assumed that the Al^O^ particles 
have been chlorinated to a depth x cm, resulting in the production 
of a hygroscopic water soluble material, the average molecular 
weight of the ions produced at complete dissociation in a dilute 
solution (M) being about 30. If a unit mass of Al20^ combines 
with surrounding gas to form a mass (1+y) of the chlorinated sol- 
uble compound, after chlorination has occurred, a particle of 

A1«0^ with an original radius of r will contain a mass p (r -x) ^ 

^ o p 3 p p 

grams of unaltered A1_0_, together with -l-Tr p (1+y) [r ^ - (r -x)^1 

dp p 

grams of the soluble product of the chlorination reaction, p^ being 
the density of the AI 2 O 2 particles. On exposure to high relative 
humidity, a haze droplet will be formed which will consist of a drop- 
let of solution enclosing a particle of AI 2 O 2 of radius . 

If this haze droplet has a radius r(cm), assuming that the 
solution is dilute, the molar ratio of the solution (y) w^ill 
be given approximately by: 


„ ^ Pp <l+y) 

I ~~Z ' moles per mole. 

M Pj- [r -(rp-x)^] 

v;here is the molecular weight of water, M is the effective 
molecular. weight of the solute (assumed to be 3a), and p, is the 
density of the solution which is taken to be equal to that of water 
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As the relative humidity increases, the haze droplet radius 
will increase. ^'Then the air is close to saturation (either slightly 
supersaturated or slightly under saturated) the equilibrium super- 
saturation of the haze droplet (S (r) ) is given approximately (in 
absolute units) by: 


S(r) 


2cr 

o 

RTpj^r 


y (r) 


S (r) is stationary (and a maximum) when 
2o M, 


RTp^r 


b _ _ ^ 
2 yr 


^ ”pPp 

3 3 ^ 

M pj^ [r -(r^-x) ] 


The 

solving 


corresponding "critical" value of r (r ) is thus found by 

c 

the equation : 

, 3 p (1+y) [r ^ - (r -x)^3 RT r^ 

f (r) = r^ - H 2 g 

2 Ma 


- 2(rp-x)^ r^ + = 0 

Since f (r) is positive at r=0, and can be shown to have one 
and only one turning point in the region r>0, f(r)“0 has either no 
real roots, or two. Since at r= (r -x) , f (r) is negative, there- 

Jr 

must be two .real roots greater than zero, one greater and one less 

than (tp-x) . Only the larger root has physical meaning, since r 

necessarily exceeds (r -x) . 

P 



The corresponding critical supersaturation (S ) is then 

o 


given by: 


2a M 

- S(r^) ~ 

L c 


In view of the definition of r , [r - (r -x)"^] is propor- 

C' c 

2 

tional to r ^ and it results that: 
o 


2a , 

c = o L_ 

c RTp^ r 2 

L c r 

c 




3 RT Mp, 


ThuS/ in a: case where the entire partic3.e has been transformed 


into a soluble material (x=r ) , it results that: 

P 

2 32 M ^ 

^ 27 Pp{l+y) p^^ RVr^^ 


which agrees with the usual formula for the critical super- 
saturation of a soluble particle v;ith a density equal to pp(l+y). 
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An evaluation of the effect of the diluted SGC in increas- 
ing the concentration of cloud condensation nuclei would there- 
fore involve solving equation (1) to find and hence evaluat- 

ting from equation (2) . A much simpler approximate procedure 
is to compare the actual SGC aerosol with one containing the same 
mass of soluble material in each particle, but without the insol- 
uble component {of radius (r^-x)) . It can be shown that the effect 
of the insoluble component in the SGC aerosol is to lower the 
critical supers atur at ion of the particle. Thus the "soluble 
component only" aerosol will provide, a lower bound to the effect 
of the SGC on the S^ spectrum and hence on cloud formation. 

Consider a haze droplet of radius r formed on an aerosol 

particle consisting of a mass m of soluble material of molecular 

weight M and containing an insoluble particle of radius r^. Then 

the molar fraction comprised by the solute is approximately 

3m Mo 

47Tp.j.j^ ' 


and hence for this composite particle,. 


2 CM 

S{r) = ° 


3m M 


RT r . ,, , 3 3 . 


(4) 


3 S 

As before S (=S(r )) is found by solving •5— = 0 to find r , 

C C d X* C 

and substituting this value of r back in the expression for S(r). 
Now, from equation (4) , 


as 

c 

ari 


, 3S(r) > . as(r) > 

^ 3r ^r 3r. '3r. ^r 

Cl 1 c 

, 35 Cr ) . 

'■a r. ■'r , 

1 c 
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since = 0 

- 3r r 

c 


Hence , 


!!c 

8r . 
c 


-9m M r . 
o 1 


1 3 2 

47rp^ M(r - r. ) 

L c 1 


Similarly, 


^ 1 / r 


3S 


3r. 9r. 

1 1 


•18m r 

i 


4ttp^ M(r 

c 


54 m M r. 

o 1 


3, 2 

> 


/i / 3 3,3 


Thus, using the fact that ( — 


is zero, it can be shov/n 


that not only the first derivative of S with respect to r is 

c o 

negative, but all the higher derivations also. Hence, the critica 
supersaturation of the particles is a monotonically decreasing 
function of r^^. 

The values for the "soluble component only" aerosol may be 

2 3 

deduced from equation (3) by substituting 3r x for r , or alter- 
’ P P 

natively by differentiating equation (4) with respect to r (v/ith 


r^ = 0), substituting this value back into (4), and finally sub- 

2 

stituting 4TfPp(l+y) r^ x for m. 

Either of these procedures yield: 

. o 32 M M 

„ 2 _ o 


81p p ^ (l+y)R^T^r ^x 

ir' J-* t' 


so that one may write 
b 


^c = r 


p 
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where b = 


32 M M 

o 

81p PT^(-l+y) R^T^x 
- p Jj 


At 20‘’C, taking y = 1, for x = 10 ^ cm, b = 4.514x10 For 


X = 10 ^ cm, b = 1.427x10 


Warm cloud effects are dis- 


cussed below for the two cases. 

The concentration (N(S )) of particles with critical 

o 

supersaturations less than is equal to that of particles 

exceeding (^) in radius. In the size range 0.035 pm<^r <25 lam, tl: 
c ^ 

distribution is given at T + 3 hours by; 

-2.5 

N(S^)= 5.938x10"^^ (|-) 

c 

while in the size range 0.01 ym<rp<0.035 ym. 


-0.75 


N(S^) = 0.2097 (^) 


The value of b depends on that of x. 

— 8 

Thus, for X = 10 cm, 

N(S ) = 1.372x10^ S 
c c 

-2 

for values of S up to 1.290x10 . For greater values of S , 

c c 

N(S ) .= 6.771x10^ s 

c c 

- 7 

On the other hand, for x = 10 cm, 

N(S ) = 2-441x10® S 

c o 


for values of up to 4.077x10 . For greater values of S^, 

ms ) = 1.606x10® S 
o o 
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The spectrum of the critical supersaturations (S } of the 

o 

natural aerosol in the Florida region has been studied by 
Fitzgerald (1972) who found that over the ocean, the average 
distribution was 

N(s„) = 2420 
^ c 

where is expressed, as earlier, in absolute units (not percent 

and N(S^) is the concentration of particles with critical super- 
saturations less than S . 

c 

Over the land, Fitzgerald found: 

N(S^) = 5913. S 
' c • c 

Depending on the nature of the aerosol and the speed of updra 

cloud droplet formation occurs on particles with critical' super- 

-3 -2 

saturations typically in the range 10 to 10 

The following Table IV-1 compares the artificial CCN concen- 
trations at T + 3 hours with the natural background concentration 
found by Fitzgerald (1972 ) . 

Table IV- 1 


Concentrations of artificial CCN activated at various supersatura 
tions characteristic of cloud formation in the SGC cloud at T + 3 
hours, and in the natural atmosphere in Florida, 


Sc 

(absolute 

units) 

-3 

Concentrations (cm ) 

-8 “7 

For x=10 cm For x=10 cm In the Natural Atmosphere 

(artificial 

particles .only) 

Over Ocean 

Over Land 

10“^ 

4.3x10° 

7.7x10^ 

1.0x10^ 

1.5x10^ 

3x10"^ 

5.8x10^ 

1.2x10^ 

1.7x10^ 

2.7x10^ 

6x10^^ 

3.8x10^ 

3.5x10^ 

2.3x10^ 

3, 9x10^ 

10"^ 

1-4x10^ 

5.1x10^ 

2.9x10^ 

5.2x10^ 
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The following Table IV-2 shows how many hours must elapse 
before the concentration of artificial CCN becomes equal to or 
less than that characteristic of the background aerosol over the 

land. When this occurs, the total concentration will be twice 

that of the natural atmosphere. The decrease of the artificial 

component has been calculated on the basis of the linear increae 

in cloud volume 'discussed in Chapter III. 

Table IV-2 

Time in hours beyond T + 3 hours until the SGC is diluted to the 
degree that the total concentration of CCN (artificial and 
natural) is twice that characteristic of the Florida peninsula. 


(Absolute 

units) 


Time (hours) beyond T + 3 hours 


X = 10“'® cm 


X = 10 ^ cm 


-3 


10 
3x10* 
6x10 
10-2 


“3 


0 

0 

2 

11 


0 

10 

24 

26 
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E, ,CONCLUSIOKS MD KESEAECH RECOMMETOATIOWS 

As .discussed in paragraph C, the concentration of giant part 
cles in the SGC is far too lov7 to result in enhanced v/arm rain 
formation. 

However, at T + 3 hours and for some hours afterv;ards , the 

SGC may contain a concentration of CCN which is large compared 

with that which is naturally present. Thus if chlorination pene- 

-7 

trates to a depth or x = 10 cm in the AI 2 O 2 particles, the SGC 

V7ill have a definite effect on cloud formation for a day or more, 

resulting in highly continental clouds which will be quite ineffi 

cient in forming rain by the v;arm rain process. In the case v;her 
"8 

X = 10 cm, this effect is probably appreciable for a few hours 
only. 

Since the effects of the SGC on cloud formation seem likely 
be very significant, it would appear to be important to make dire 
measurements of the spectrum of the critical supersaturations of 
SGC aerosol. This would provide in a direct manner a far more ce 
tain basis for the estimation of warm cloud effects than is pre- 
sently available, or could be deduced from chemical analyses. 
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Chapter V 

CLOUD PHYSICS PKOCESSES - COLD CLOUDS 

A. ICE PHASE DEVELOPMEM' AND PEECIPITATIOH IN CLOUDS 

1. General 

Only a small percentage of clouds reach the precipitation stage. 

The progression from minute cloud droplets (circa 1-25 pm radius) to 
falling hydrometeors 3 if it is to occur at all, involves three basic 
processes : 

a. droplet collisions and coalescence 

b. the ice crystal or Bergeron-Eindeison process whereby ice 
crystals grow by diffusion of water vapor at the expense of 
evaporating supercooled drops and from cloud vapor generated 
in vertical updrafts 

c. the ice crystal process augmented by collisions with droplets 
(riming) and/or other crystals (aggregation) 

Warm rain or that due entirely to droplet coalescence dominates 
at tropical latitudes. It may even play a role at higher latitudes with 
unstable clouds not extending to or far above the freezing level. The 
ice crystal mechanism clearly is dominant at polar latitudes and also 
highly significant at mid-latitudes. In the latter zone, where the 
world's population and industrialized nations are concentrated, the 
ice phase combined with collisional mechanisms (item c) prevail. It 
is well recognized that the heavier mid-latitude precipitation (rain or 
snow) can only be explained by this combination of mechanisms 
( Houston , 1950 ) . 

2. Florida Precipitation Processes 

Florida, being a sub-tropical region, can experience rainfall by 
either mechanism (a) or (c) above. However, deep cloud systems and the 
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ice phase undoubtedly are instrumental in the major production of rain- 
fall on the peninsula. As described in Chapter VII, the sxunmer rainy 
season extends from roughly May to September or October. During this 
time, rainfall is likely every day (505^ probability) vs. 1-2 days per 
■week in winter, and half the rainfall comes from local shoi>rers and 
thundersho'wers (Bradley, 1972). 

Clearly in these deep convection systems , ice nuclei and crystals 
are the initial building blocks for subsequent riming, snowflake 
aggregation, latent heat release, and hea-vy rainfall. The Florida Area 
Cumulus Experiment (FACE), conducted by HOAA since 1970, is predicated 
on the belief that cloud seeding with additional ice nuclei in organized 
convective systems •will enhance rainfall. The effectiveness of FACE to 
date is not entirely clear; however, the recognized role of the ice 
phase in Florida precipitation and the comparison between ice nuclei 
seeding concentrations and that inadvertently released in NASA, space 
shuttle launches are highly relevant (Section C). 

Standard temperature lapse rates for summer and winter at latitude 
30°IT (AFCEL Handbook of Geophysics, 19^5) are sho™ in Figure V-1. 

It is e'vident that the freezing levels during the respective warm and 
cold seasons average roughly 5 and 3.5 Im. The -5°C levels, below which 
ice crystals are only rarely observed, are at approximately 6 and 4 km. 

As a first approximation then, cloud systems with tops below these levels 
will not. be influenced by ice and only warm rain processes will be involved. 
In strong convection, particularly in summer, Florida clouds can pene- 
trate well above 6-8 km where freezing phenomena are paramount. 
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Relative Humidity (%) 


Temperature ('“K) 


Figure V-1 Subtropical Lapse Rates (30"N) 
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B. ICE MJCLSATION MODES MD EATUBAL ICE EUCLEI (IE) COECEMRATIOES 

1. Ice Eucleation Modes 

Particles that promote. the formation of ice in clouds are believed 
to do so via: 

a. condensation-freezing 

b. contact (with a supercooled drop) 

c. sublimation (or direct deposition of vapor onto a solid 
nucleant ) 

Ehile imcertainty remains , the condensation- freezing process is a 
principal mode of nucleation in the atmosphere. In short, mixed ice 
nuclei, consisting of mainly hydrophobic composition with some hygro- 
scopic sites are effective in attracting water and then initiating 
freezing. Evidence suggests that this nucleus type may represent the 
most general type of IE. Contact nuclei, necessarily very small and 
hydrophobic to avoid building up a water film, appear capable of 
freezing contacting droplets at relatively warm temperatures. Sublima- 
tion nuclei now are considered to be comparatively rare in natural cloud 
processes, although artificial seeding agents such as silver iodide can 
act in this manner. 

2. Eatural. IE Concentrations and Sources 

'^'Jhile condensation nuclei are plentiful in the atmosphere, parti- 
cularly over continents and at mid-latitudes, ice nuclei are scarce. 

The well-quoted global concentration of Ifl--is lt“^ at -20°C (some 6 
orders of magnitude less plentiful than cloud condensation nuclei ) . 

For every U°C temperatiore warming, there is an approximate order of 
magnitude decline in activated ice nuclei. Figure V-2 illustrates the 
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average concentration of lU as a function of temperature as obtained 
from the en^iirical function of Fletcher (1962 ): 

IH (A“^) = 10-^ exp (0.6 at) 

where AT is the degree of supercooling. While order-of-magnitude 
departures from this cu 2 rve occur, it typifies measurements of natural 
m concentration made to date. 

During the last 10-15 years, however, observational evidence has 
accimulated to the effect that ice crystal concentrations in clouds 
sometimes greatly exceed that expected on the basis of ice nuclei con- 
centrations — by factors of up to 10^ . Such occurrences seem more common 
in aged maritime clouds containing some relatively large drops (e.g. , 
Mossop, 1971). Recent laboratory results (Hallett and Moss op, 197^) 
indicate that this so-called ice-multiplication process may be at least 
partially resolved. Crystal enhancement apparently necessitates riming 
of drops larger than 25 ym dia. at selective temperatiures very close to 
-5°C. In a summary of ice crystal (IC) measurements in clouds by 
several investigators, it was generalized (Hobbs, 197^) that the IC/IH 
ratio was on the order of 10**^ at decreasing to about 10^ at 

-lii°C and 10^ at -20°C. 

Thus in those selective clouds where multiplication has taken place, 
background levels of ice might be of order 1 at — lU°C and 10—20 
at -20°C (in- contrast to Figure V-2 values). The foregoing meas^lrements 
and estimates will be of relevance in estimating the importance of 
potential IH released in solid rocket motor (SRM) launches by HASA 
where maritime clouds are often prevalent. 
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While the sources of naturally occurring ice nuclei are also tenuous, 
the activity spectra of several suspected materials have been analyzed. 

The earth’s surface, a logical source, contains clays, silicates, and 
minerals that can serve as active IW (Schaefer, 19^9 i Mason and Mayhank, 
1958 ; Mason, 1971 ). Clays have activity thresholds (l active nucleus in 
^'10‘^) at temperatures as warm as -10°C (common kaolinite at -9°C) and 
reach high activity levels hy -24°C. A list of some of the materials 
examined and their ice nucleation thresholds appear in Table V-I. Prom 
collected snow crystals which were then sublimed, Kumai ( 1961 ) observed 
that soil particles (clays) were the commonly found central residue. Of 
relevance and to be discussed subsequently, the stabilized ground cloud 
(SGC) contains large amounts of soil material sucked up from the ground 


d^lring rocket laimch. 



Table V-P 
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Ice Nucleating Ability of Naturally Occuring Particles 
(from Mason 1971) 


repeoducibility 

AT T> A AW TR POOE 


Substance 

Chemical 

cempesition 

Symmetry 

Threshold 

temperacuro 

("0 

Boforencos 

CoVoUito 

CiiS 

Sex 

—5 

.M*M 

Vatcrito- 

CaCOa 

Eox 

-7 


^•Tridymite 

SiO, 

Hex 

— 7 

JI&H 

Magnetite 

Pe^O, 

Cubic 

-s 

M&il 




-D 

I&I 

H&otiuite 


P^ncUolc 

-9 

M&ai 




-13 

i&i 

Annuxite 


Monoclinic 

-9 

31 

lUito 


Monoclmic 

-9 

31 

Metabontomte 



-9 

31 

Microclmo 



-9 

31 i3I 

Hyporsthene 

(Mg,Fe).(S<,Os) 

Hhomblc 

-10 

I&I 

Hoomutito 

Fo.Oj 

Hex 

-10 

3I&3I 

(Spocularite) 



-13 

I&I 

PyTTophyllita 


Monoclinic 

-10 

31 

Gibbstto 

A[;.{0H)5 


-11 

M Ss2,l 

HaUoysito 

AJ,(0IJ)4 Sirt- 





Oj-2S,0 

Monociinic 

— 13 

M 




-13 

3I&3I 

Du^kitc 

Al»(OH)., 5i.O, 

Jlonodinio 

-13 

31 

Olivine 

(Ug, Fa),-S!0^ 

Hbombic 

—12 

P &S 




-IS 

I&I 

Aquadag 

c 


—13 

M 

Dolomite 

CaMgfCOj), 

Hex (Ehoinb) 

-14 

31 &3I 

Biotito 


Monoclimo 

-14 

31 &31 

Attapulgjto 

4H.O-(OE)„Mgs- 






Monoclinio 

-14 

31 

Muscovite 


Monoclmie 

-14 

I &I 




<-is 

31&3t 

Vermiculjte 


^ionoclimo 

-la 

M &3I 

Xoatroaite 


Monochnic 

—15 

31 

MontznoriUonite 


Monoclmic 

-16 

31 




<-IS 

31 &31 

Gypsum 

CaSOi'SHnO 

llonooliais 

-16 

3I&3I 

Cranbite 

0 


-16 

31&3I 

Cixmabar 

HgS 

Hex 

-16 

3I&31 

OrthDcIesa 

KAlSijOs 

Monociinic 

—17 

IiKl 




<-18 

M&3I 

Anortboclaso 



-17 

31 &31 

Quartz 

SiO, 

Hex 

<— IS 

3t&3t 




-19 

I&I 




A 

.1 

O 

r,K&o 

Stony metcorito 



—17 

3I&31 




-IT 

I&I 



(2 speennens) 

<-is 

3I&31 



(3 Specimens) 

<-17 

yi 

Volcanic ash: 





10 Japanese volcanoes 


-12 to —16 

1.K&0 

JIt. Etna 



-13 

31&3I 

Crater Lake, Oregon 



-10' 

S 

Parienten, ^loxico 



-33 

S 

Soils: 





Arnerican-lcom, clay, 


—3 to —35 

S 

clay 



— 11 

P &S 

Loess, N, China 



-12 

I&I 




— 15 

I,K&0 

Looss, IInnford,^ir.SA. 


-11 

S 


Tho following substances voro found by Mason d al. to bo macUvo at tomperatnros 
above — 18 ^ 0 ; Montmorsllonite, ssptoHtOj albltOj talc, sand, quartz, c(‘ind>*Tnite, aad 
sovera! samples of stony moteortto. 

M (StM ^lason and Maybanlc (195S), M = Mason (l9C0c), 1^1:= Xsono and Ikebo 
(1060), X, K & 0 Isoiio, IComabayaii, and Ono (1950),, P & S = Truppacher and 
Sanger (1955), S = Schaefor (19i95). 
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C, ICE HUCLEATION OE FLOEIDA CLOUDS BY THE S.G.C. 

!♦ ActiTity of Rocket Exiia\ist Particles - AI 2 O 2 
As noted in Chapter III, the stabilized groimd cloud is generally 
confined to the first 2 km above ground. Within that altitude increment, 
the corresponding rocket burn yields approximately 7-*10 x 10^ g of 
see sample Table V-II (from NASA JPL Tech. Memo 33-712). 

Some uncertainties exist as to the dominant crystalline phases and 
exact surface composition of the AlgO^ particles, particularly upon 
interaction with gaseous and aqueous HCl. Investigators at NASA-Langley 
Research Center (LaRC) (Gofer and Pellett) have suggested that a major 
proportion (by number and also surface area) of the alumina produced by an 
SRM tends to consist of submicron-size metastable alumina types. This 
judgment was based partly on LaRC x-ray diffraction resvlts obtained on 
alumina collected by JPL from small SRM's fired in a closed tank; findings 
agreed with the electron diffraction results of Robbins and Strand (l970), 
i.e., the large spherical peirticles (>0.2 to 1 inn) were predominantly stable 
alpha- alumina and the irregular submicron particles (highly agglomerated 
as received) were predominately metastable gamma-alumina. 

However, the slower effective particle cooling rate in the exhaust 
plume of a large SRM may significantly increase the percentage of the alpha 
form, as observed in recent measurements by Varsi (1976). He reported that 
Titan III particle samples , collected at high altitude by NASA-Ames and 
examined by both transmission and reflection electron diffraction at JPL, 
indicated "...mostly alpha phase for particles in the 0.1 ym size range." 
Although earlier JPL transmission electron diffraction measurements 


"...on just a few 0.1 ym particles collected by JPL with an impactor" at 



V-10 


Table V-II Exhaust Products Released into Atmosphere 

by Space Shuttle SRMs (Partial list: Mission 2)* 


Altitude band, 
km 

A a Iti tilde , Time, 

km s 

Atime, 

s 

Average 
mass 
flew, 
10^ 5<g/s 

Mass, 

1 0^ kg 


Sxhaust products^ 



AI2O3 

0. 3020 28 

HCl 

0. 209315 

CO 

0.241719 

Exit 

CO, 

0.0343‘J46 

a.xit 

0-0. 0095 

0. 0095 

0-2 

4 

> 

9.444 

18.88 

5.704 

3.953 

4. 565 

0. t>495 

0.009S-0. 039 

1 

0. 0295 

2-4 

i 


9.446 

18.89 

5.704 

3.953 

4.565 

0. 649b 

0.039-0. 087 

0. 048 

4-6 



9.446 

i 

5.705 

3.954 

4.566 

0.6496 

0. 087-0. 160 

1 

0.073 

6-8 

i 


9.447 

■ T 

5.705 

3.954 

4.566 ■ 

0.6497 

0. 16-0.25 

0. 090 

8-10 

2 

9.448 

18.89 

5.706 

3.953 

4.566 

0.6498 

0.25-0. 50 

0. 250 

10-14 


1 

9.449 

37.79 

11.41 

7. 912 

9. 134 

I.; 

500 

0. 50-0. 85 

0. 35 

14-18 



9.451 

37.79 

11.41 

7.911 

9. 136 



0.85-1.3 

0. 45 

18-22 



9.451 

37. 80 

11.42 

7.912 

9.137 



1.3-1. 9 

( 

0, 60 

22-26 

4 

9.454 

37.81 

11.42 

7.914 

9.139 

1. 300 

1.9-2. 2 

0. 30 

26-28 

2 

9.445 

18.89 

S.704 

3.953 

4.565 

0.6490 

.2.2-2. 5 

0. 30 

28-30 

2 

9.292 

18.58 

5.612 

3.889 

4.491 

0.6390 

2.5-3. 3 

( 

0. 80 

30-34 

/ 

1 

8.859 

35.43 

10,70 

7.416 

8.564 

1.219 

3. 3-4. 2 

( 

0. 90 

34-38 



8.282 

33. 12 

10. 00 

6.933 

8.006 

1. 139 

4.2-5. 1 

0. 90 

38-42 


1 

7.705 

30.81 

9.307 

6. 450 

7.448 

1.060 

5.1-6. 0 

0. 90 

42-45.44 

- 

1.44 

7. 123 

24.50 

7.399 

5. 128 

5.921 - 

0. 8426 

6-9 



45.44-56.06 

10.62 

6.526 

69.29 

20.93 

14.50 

16.75 

2.383 

9-12 



56.06-65.52 

9.46 

6.630 

62.70 

18.94 

13. 12 

15. 16 

2. 157 

12-15 - 



65.52-74.01 

8.49 

6.925 

58. 78 

17.75 

12.30 

14.21 

2.021 

15-18 



74.01-81.64 


’.63 

7. 147 

54.52 

16.47 

11.41 

13.18 

1.875 

18-21 



81,64-88.70 


’.06 

7.270 

51.32 

15. 50 

10.74 

12.40 

1. 765 

21-24 



88. 70-94. 90 

6 

1.20 

7.320 

45.37 

13.70 

9.497 

10.97 

1. 561 

24-27 



94.90-100,8 

5. 90 

7.298 

43.05 

13. 00 

9.011 

10.41 

1.481 

27-30 



100,8 -106.2 

5 

.. 40 

7-227 

39.02 

11.78 

3. 167 

9.431 

1. 342 

30-33 



106,2 -111.2 

5. 00 

7.039 

35.44 

10.70 

7.417 

8.565 

1. 219 

33-36 



111.2 -116.08 

4.88 

6.230 

30.40 

9. 180 

6.362 

7. 347 

1.045 

36-39 

3 

116.08-120. 80 

4.72 

3,227 

15.24 

4.602 

3. 189 

3.680 

0. 5241 

39-41.6 

2,6 

120. 80-124.85 

4. 05 

0.8686 

3,519 

1.063 

0.7366 

0.8506 

0. 1210 





TOTAL 

915.6 

276.5 

191.6 

221.3 

31.49 

The mass is indicated below the symbol for each specie. 


REPRODUClBILrH OF® 

ORIGINAL PAGS H 


■'NASA JPL Technical Memorandum 33-71 
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lower altitude, showed the gamma (cuhic) stmcture, Varsi (1976) concluded 
that "At present it would appear that the evidence_ is in favor of alpha 
(hexagonal) structure for the material injected into the stratosphere." 

A recent assessment of the interaction of various aluminas with gaseous 
and aqueous HCl (Cofer and Pellett) indicated that both the metastahle (theta, 
delta, gamma) and stable alpha forms chemisorb gaseous HCl, either dry or 
moist, to yield significant coverage of the surface by soluble chloride 
("saturation" at 12 A^molecule) . In fact, recent studies of anhydrous HCl 
adsorption (up to 20 torr) at Virginia Polytechnic Institute (NASA-LaRC Grant) 
have indicated that surface coverage for an alpha sample (olT A^/moleeule) 
was about twice that for a Alon-C gamma sample ('\j 33 A^/molecule) . 

Summarizing possible AlgO^-HCl interactions , while the metastable aluminas 
tend to be more soluble in the bulk than the alpha (hexagonal) variety, all 
the known forms of alumina appear to chemisorb HCl and exhibit the potential 
for attaining relatively high surface coverages of soluble chloride. Thus 
the probable result of in-cloud chemisorption processes is the production of 
partially soluble particles which tend to be less effective as ice nucleants; 
how much less is uncertain. 

Since alpha AI 2 O 2 is a reasonably active lU source and because of phase 
and surface-composition uncertainties , we will tentatively assiome that as 
much as 1 to 10^ of the total AI 2 O 2 is of the stable type. As subsequent 
estimates will show, the exact percentage is highly important and must be 
determined acciarately. 

Thanks to the cooperation of the China Lake Haval Weapons Center and 
Drs. William Finnegan and Kirk Odencrantz, nucleation tests were run at their 
facility. They burned an aluminum solid-propellant in their 15 m^ cold chamber 



V-12 


and determned the temperatiare a’ctivity spectrum of the stable so 

produced. Their results appear below in Table V-III. 


Table V-III. AI2O2 (Stable) Ice Uuclei Activity 


Temp. 

-lU to -15°C 


-20 


Output per gram 
1-2 X 10^ g“l 
A.1 x,10l° g-1 


(Agl comparison) 
10^^^ - lO^S g-^ 

10^5 - 3 X 10^5 g-^ 


In short the threshold of activation was at -l4 to -15°C with the 

IN concentration reaching '^'1 x 10^® g“^ at -20°C. For reference, the 

typical effectiveness of the best cloud seeding agent, Agl, is also 

indicated for various generators (Fletcher, I962). 

The concentration of effective HJ in the SGC is at any time (cloud 

volume) and temperature given simply by; 

Mass of a 1202 x Activity (gT^) 

*" Cloud Volume 

For example, at (T + 3) hours the cloud volxmie is 6 x 10^ 

(6 X 10^*^ t). Assming the unreasonable case where all the AlgO^ is 
nonchlorided (stable) and T = -li^°C, then: 

lU = (7 X 1Q7 g)(lQ8 g-1) 

6 X 10^*^ £ 

A.djusting for the perhaps still conservative values of 1—10^ stable 
AI2O25 the concentration of IN at -l4°C woxxld be 0.1 - 1 £“^. Continuing 
in this manner for two cloud temperatures ; the 1-10^ stable AI2O2 
assumption; the expanding cloud volume with time and initial AI2O2 mass 
of 150 yg as given in Chapter III; the IK concentration values of 


Table V-IV result 
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TaMe V-I¥. Ice nuclei Concentrations from AI 2 O 3 (l-lO^ Active) 


(Values versus time and temperature) 


Time 

(T + 3 hrs ) 
(T + 1 day) 
(T + 3 days) 


Cloud 

Voluae 

6 X 10^ km^ 

5 X 10 ^ 

1.5 X 10*^. 


■ IJT (~l4‘=*C) 
0.15 - 1.5 


.02 - 0.2 


.006 - .06 


in (-20°C) 

15 - 150 
2-20 
0.6 - 6 


CoD^aring Table V-IV estimates of potential AI 2 O 3 ice nuclei with, 
the natural atmospheric m background of Fig. V-2, it is evident that the 
added component is not insignificant. If 10J^ of the AI 2 O 2 is of the 
stable type, the ice nuclei conceivably produced is roughly 1-2 orders of 
magnitude above background at (T + 3) hours; and still a factor of 2-5 
above background after 1 day. Even for stable AI 2 O 3 and (T + 3) hours, 
note at -20°C that 15 is well above the often quoted natural atmos- 
pheric level of 1 In maritime clouds where ice multiplication may 

be tahing place, the differences between SGC and natural IIT concentrations 
are not so great. After three days natioral IE levels generally (but not 
necessarily always) exceed the values in the expanded rocket cloud volume. 
Precipitation scavenging, during the rainy season in particular, would 
likely reduce the residence time and concentration of SGC particulates 
to insignificant values at ''^3 days and longer. It must be emphasized 
that this entire rationale hinges on 1 - 10 ^ of the AI 2 O 2 being effective 
as IE and the rest not at all. 

2. Ground Component of Ice Nuclei 

An added consideration is the amount of ground soil drawn up into 
the SGC dtiring rocket launch. If Fig. V-3 (NASA Notes Ref. l4) for a 



Concentration of, Suspended Particles, /ig/m' 


1000 


100 


lo' 


10 100 500 

Downwind Distance from Launch Pad, km . 


Particles in Titan III Rocket Cloud Sampled at Ground Level 
(HASA Notes, Ref. l4) 


1 


Figure V-3'‘ 
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Titan III launch is representative, then the total suspended material sampled 
at ground level is some 30-100 times that of the AI2O2. Some of this ground 
zaaterial will undoubtedly constitute potential ice nuclei. Most clay- 
silicates, as shown in Table V-I, will be more effective (warmer nucleation 
thresholds) than the AI2O2, provided no HCl "poisoning" takes place. These 
active particulates are often found in trace amounts in the eaxth’s crust 

and may not be abundant at Cape Kennedy; typical sand is not active at 
temperatures above -l8*^C (Mason, 1971) • Nevertheless, it is inevitable 
that some ground particles will serve as IN and the concentration may 
well be significant, as the previous numbers demonstrate. More research 
and instrumentation on the chemistry of SGC particulates and IN 
activation spectra are clearly in order. 

3. Cloud Seeding Implications 

Thus the potential for some cloud seeding exists, particularly 
within several hours or perhaps a day of launch. The higher values of 
Table V-IV (l5 - 150 5-”^) represent substantial seeding concentrations; 
cloud seeding programs often are designed on comparable values. The 
NOAA Florida FACE program, aimed at "massive" cloud seeding, attempts 
to inject 100 Agl nuclei active at -10°C levels. An example of the 
effectiveness of small piyotechnic seeding rockets (U.S. Navy and Olin 
Corp. ) used on the Florida program is shown in Fig. V-U. While their 
Agl activity (g”^) at -lii to -20°C is 3-4 orders of magnitude higher 
than stable AI2O3, the great quantity of AI2O3 released in shuttle 
rocket burn results in the relatively high IN concentration of Table V-X7. 

Returning to the impact on Florida peninsula clouds in summer, the 
AI2O3 activation threshold of -l4°C corresponds to an average 7.2 km 




Figvire-.y-U Nucleating Effectiveness of Seeding Rockets (NOAA. FACE Program; 
Woodley and Sax, 1976) 
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aQ.titucLe (Fig. V-l). Note that temperature structure can depart sub- 
stantially from the standard lapse rate. Hoirever, in general it can 
be stated that clouds vith tops below 6-7 km •will not be influenced 
by AI 2 O 2 ice nuclei (in winter, 4. 5-5. 5 kn cloud tops appear tolerable). 

In deep penetrative convection associated with well developed cumuli 
and thunderstorms , seeding material most assuredly can reach levels of 
substantial activation (roughly 7~10 km in summer j 5-9 km in winter). 

Above these levels (-35°C) , generally there are natural IN concentra- 
tions far greater than man can produce. In fact at approximately -UO°C 
(lO-ll km), spontaneous crystallization of supercooled water drops 
occurs ■without the need of any nucleating particles. 

The significance of enhanced IK concentrations is far more diffi- 
cult to assess. Exact cloud seeding effectiveness has been and continues 
to be a subject of debate. In broad terms, given suitable environmental 
conditions and substantial supercooled clouds , IN seeding of the order 
of 10 SL~^ is believed by some to increase precipitation by perhaps 
10-20^. More massive seeding (circa 100 fi.“^) in thunderstorm airmasses 
reportedly can diminish damaging hailfall (Sulakvelidze, et al. , 1967; 
Burtsev, et al. , 1973; Miller, et al., 197^). These t-wo weather modifi- 
cation effects — potential rainmaking and/or thunderstorm diminution — are 
most relevant to Florida. Neither are necessarily detrimental, especially 
the latter. Alternately, seeding at an inopportune time can have the 
effect of suppressing cloud development and rainfall (Braham, I966). It 
has been postulated (Schleusener, I968) that Agl seeding rates higher 
than 2 kg hr“^ tend to enhance convective activity. 
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k. Eocket Release Above the SGC 

While this analysis is concerned, with the low-level SGC, the 
possible IWM effects at higher levels shoiild not be ignored. It is 
evident from Table V-II that approximately 8 x 10 ^ g of AI2O3 are 
released at 2-12 km altitudes, a mass eq.ui valent to that in the first 
two km. Considering that the vertical extent is 5 times as great and 
that dispersing winds are greater at higher levels , the resultant ice 
nuclei concentrations of Table V-IV should be down by roughly an order 
of magnitude at these higher altitudes. Also ground soils, ingested into 
the SGC, will not immediately reach these levels except by subsequent 
strong convection. Thus, the possible IWM ingiact by direct AI2O3 
release at high levels, while a potential factor in certain circum- 
stances, is proportionately diminished. 

5. Summary IWM Implications 

In summary, on the assumption that 1 - 1 . 0 % of the space shuttle 
rocket AI2O3 (and/or entrained earth material in the SGC) are effective 
ice nuclei with a threshold of -lij-°C: 

a. The potential for inadvertent weather modification (IWM) exists. 

b. The effect coTild be that of altering precipitation amounts, 
hail, and severe winds; in the uncontrolled situation involved, 
the net result could be either an increase or decrease. 

c. Concerning rainfall it is more likely that such an effect would 
lead to an increase of modest amount and be of modest signi- 
ficance (based on non-orographic cloud seeding conducted to 
date); because of the crucial timing and sizable seeding re- 
quired to modify hail development , significant alteration 
appears more improbable, though possible by chance. 

d. Seeding effects are more likely in summer when strong convection 
can carry particulates upward to colder -IN activation levels. 
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e. The levels most conducive to ice nuclei crystallization are 
approximately 5-10 km, the higher end of the range in summer 
and the lower levels in winter. 

f . Any rWM is more probable at shorter times (T + 3 hours ) , owing 
to higher UJ concentrations, with the impact diminishing with 
time. Concentrations may still be somewhat above background 
after one day in continental type clouds but probably not 
enough so to perturb weather significantly. At 3 days and be- 
yond, IWM is considered highly improbable. 

g. released above the SGC in the 2-12 km altitude range 
are less concentrated by about an order of magnitude. Some 
near-term short-range IVM could result if susceptible clouds 
are present. 

h. Because of washout and dilution effects of the SGC with time 
(particle residence time of a few days in the lower tropos- 
phere), no cumulative Xl^M effect from the projected liO launches 
per year is likely. As an added precaution, spacing of rocket 
launches by several days is recommended. 
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D. BISK ANALYSIS 

1. Evaluation Guides 

The previous section (C-5) referred specifically to the weather 
modification risks associated with shuttle rocket launches. In parallel 
fashion, one can examine those conditions most favorable to the WOAA 
FACE seeding program. Their suitable criteria for rain augmentation 
in Florida are listed in the first column of Table V-Y. All ''s^^itable'' 
conditions are relevant risk guides -for NA.SA purposes ■srith the exception 
of item 7. Here, even isolated (large) cumuli could be substantially 
stimulated on a local scale. 

To this list of most susceptible seeding conditions influencing 
Florida, one should add the following; 

a. short time periods (hours) after launch 

b. when thunderstorms or large cyclonic systems are in the vicinity 
at time of rocket launch 

c. when winds are calm or easterly 

d. items (c) through (g) of the previous section relating to 
most conducive season (summer); altitudes (5-10 km); strong 
convection; and cumulative consequences (only with short- 
intervals between launches). 

Conversely, the second column of Table V-V indicates conditions least 
favorable for weather modification — planned or inadvertent . Wherever 
possible, NASA laimch schedules should attempt to take advantage of 
these mesoscale and cloud-microphysical factors. Prevailing westerly 
winds would also assure that any possible effects take place over 


the ocean and not over Florida. 
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Florida NOAA FACE Seeding Criteria (Woodley and Sax, 1976) 


Condlclons Suxcablc for 
''Dynamic” Sectllng of Florida Cumuli 

1. unscable lower 'troposphere, stable middle 
troposphere and unstable upper troposphere 
with moisture values not too dry or coo 
wet throughout 


2. hard cauliflower appearance to clouds 


3. tower moving upward through flight 
altitude 

4. cloud liquid water content in excess of 

1.0 as measured by Johnson-Williama 

device 


5. icG particle content less then 7.5 per liter 


6. updraft velocity in excess of 7.5 m aec"^ 

7. convection neither very isolated nor very 
disturbed 


8. weak wind shear conditions 


9. weak low-level winds 


10, continental air mass characteristics 


Conditions not Suitable for 
"Dynamic” Seedln;?, of Florida Cumuli 

stable lower tropoephero or unstable troposphere 
through mid- levels; very dry or very wet through 
large region of troposphere 

fuzzy appearance to clouds 

no upward motion to tower 

cloud water content less than 0.5 as tf 

measured by Johnson-Willlams device 

ice particle content greater than 20 per liter 

updraft velocity less than 5.0 n sec*^^ 
suppressed conditions or disturbed conditions 

strong wind shear conditions 

strong low-level winds 
maritime air mass characteristics 


Explanation 

For optimum dynamic seeding effects, cumulus 
clouds need to grow naturally to about -ID^C 
but should be stopped by a capping Inversion; 
heat released by seeding allows cloud to grow 
beyond inversion into unstable layer above. 

Hard appearance is Indicative of young, . 
vigorous cloud with a copious quantity of 
supercooled water. 

Cloud should be in growing stage of Its life 
cycle for optimum seeding effect. 

Johnson-Uilliams instrument measures water 
content of drops <30pm radius; the presence of 
supercooled cloud drops is indicative of a 
youthful tower which has not yet started on 
the "downhill” part of its life cycle. 

Seeding Is pointless if large concentrations 
of natural] y-formod ice already exist In 
\ cloud updraft region. 

Strong updraft is Indicative of youthful 
cloud. 

Seeding effect becomes lndistingul&/hable If 
all clouds grow to greet heights, however, 
cumulus cowers should form in close proximity 
If seeding is to cause mergers. 

\ Cirrus blowoff can decrease surface heating 
and thus suppress new convection; also natural 
seeding from cirrus anvils can confuac experi- 
mentacSiOn. 

Continuity of thermals destroyed by strong 
low-level winds and convection is suppressed. 

It is suspected (though not confirmed) chat 
natural glaciation occurs more readily under 
maritime air mass conditions. 
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2, Magnitude of Possible XWM 

There is reasonable agreement that the seeding of layered super- 
cooled clouds over orographic terrain can produce precipitation increases 
of the order of 10-20^ (National Academy of Sciences, 19T3). The 
seeding of cumulus clouds brings into play more cloud and environmental 
variables and less predictable results (Neiburger, 1969-). Perhaps the 
most positive such experiment conducted to date is the Israeli I 96 I-I 966 
program in which an average l8% precipitation increase was reported. 

Simpson and Woodley (l9Tl) reported cases of explosive growth of seeded 
maritime clouds with perhaps a doubling of precipitation at cloud base 
level as revealed by radar. (Ground level precipitation enhancement 
would logically be substantially less, but no corresponding measure- 
ments were obtained. ) There are numerous other cumulus programs that 
show little, if any, increase in precipitation and in some cases decreases. 

Based on these deductive considerations, the upper limit risk factor 
for IWM on the NASA Space Shuttle Program could be 20^ on a local scale 
and over short time periods (hours). This is considered a pessimistic 
view and except for rare situations, the risk factor should be much less. 
Should the percentage of active AI 2 O 3 in the total exhaust be less 
than 1% and should ground soils in the "vicinity of Cape ^Cennedy be 
generally devoid of effective clay-silicates, effects via cloud 
glaciation reduce to the noise level. 

Until more information is obtained on the exact chemistry and 
ice activation spectra of rocket and soil material within the SGC, 
caution is advised. Launches, whose rocket plumes could interact 
with existing or predicted large cumuli (or a strong sea breeze regime) 

i 
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•within 1-3 hours, are not recommended. 

3. Recommendations 

It is ad'vlsed that further MASA efforts he conducted to provide 
information ■'ri.tal in validating certain ice phase assumptions necessary- 
in this evaluation. These recommendations are as follo'ws: 

a. determine the percentage of stable AI 2 O 2 generated in the 
shuttle exhaust plume 

b. perform additional ice nuclei activation spectra of AI 2 O 3 
(as done by the China Lake Kaval Weapons Center for this 
analysis) and, if possible, 

1 . at high temperatures approaching rocket burn and 

2. with the identical SEM propellent mix 

c. in future penetrations of the stabilized ground cloud, include 
membrane filter collections of particulates . These filters 
can subseq.uently be processed in the laboratory for perhaps 
the most definitive indication of ice nuclei active at • ' 
various temperatures (also fly an automatic IW counter) 

d. establish a groimd network to evaluate possible downwind 
changes in precipitation patterns and to collect rain water 
for chemical analysis (e.g., AI 2 O 33 pH). 
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Chapter VI 

SOLim ATTENUATION MODEL FOR THE STABILIZED GROUND CLOUD 

The effects of aerosols on the atmospheric transmission of solar 
energy have heen debated since early in this century. Humphreys ( 1913 ) 
concluded that volcanic dust caused a cooling in the average surface 
temperature of the earth. The effect of the Increasing haziness around 
cities j which is attributed to increased air pollution, has been studied 
by numerous authors. Bryson (1968) has "concluded the steady decline in 
average temperature of the earth over the past decade is a result of 
increasing global pollution levels, Mitchell (l 9 Tl)j also working on the 
problem, showed that increasing pollution may lead to a warming or a 
cooling of the earth depending on the aerosols ' radiative properties . 

On a smaller scale, the effects of increasing air pollution in urban 
areas have been summed up by Landsberg (1962) and Peterson (1969). Their 
findings indicate that urban areas receive an average of 15 to 20 ^ less 
solar radiation than rural regions. The reports quoted here indicate the 
uncertainty in ascertaining the direct effects of air pollution on solar 
energy. 

The model was designed for rapid operation by reducing computational 
demand wherever possible. In addition, the model was designed to require 
only easily accessible input data. This criteria served the final goal 
of maintaining generality and versatility in the models when possible. 

The transmission of solar radiation through the atmosphere is com- 
puted using an iterative techniq.ue to solve the radiative transfer 

equation (Dave and Gazday). This technique first approximates the 
total intensity as that due to first-order scattering only. Successive 
approximations include the next highest order or scattering until the 
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desired acctiracy is obtained. Thus , the final calculated intensity in- 
cludes multiple-scattering effects to the order of accuracy desired. 

The transmission of sunlight through the atmosphere is a complex 
phenomenon, Eough estimates of the loss of radiation from a direct beam 
of sunlight can be obtained q.uite easily using single-scattering techniques 
(Van De Hulst, 195T). HoweTer, since ve are interested in determining a 
heat budget at the ground, we are required to use more complex techniques 
referred to as multiple-scattering approximations. Single-scattering 
techniques yield the amount of solar energy lost from the direct beam of 
smlight (direct radiation) ; the multiple scattering method, however, 
yields the single-scattering result with the addition of the light 
(diffuse radiation) scattered in the downward direction. Since the earth's 
surface collects all sunlight directed downward, the multiple-scattering 
result is significantly more accurate for our purposes (Hammond, 1973) . 

Multiple-scattering calculations are considerably more complex 
because they require the solving of an integro-differential equation 
called the radiative transfer equation (R.T.E.). The equation has 
historically been used largely in the fields of astronomy and physics. 
However, with the realization that anthropogenic air pollution may effect 
the global radiation balance , and hence the climate of this planet 
(Charlson and Pilet, 1969), the R.T.E. recently has been used to calculate 
the effects of pollution of the atmospheric radiation balance of the 
earth . 

In principle we know how to solve the R.T.E. using one of the various 
numerical techniques. The most widely used methods include the doubling 
method (Hansen, 1971), discrete ordinate method (Chandrasekhar, i960), 
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spherical harmonic method (Bergstrom and Viskanta, 1972), and the iterative 
method (Braslau and Dave, 1972). The computational demand of the dis- 
crete ordinate method is giiite large, thus this method ■was not utilized 
for the model. Of the remaining three techniques , the itera'cive method 
was chosen for this model because program modifications are easily per- 
formed using this method. 

The primary difficulty encountered in applying these techniques is 
the calculation of the scattering and absorption properties of the medium. 
Scattering and absorption of electromagnetic radiation (light) by par- 
ticles or droplets are described by Mie scattering theory. Mie scattering 
theory is the classical solution describing the interaction of electro- 
magnetic waves with spherical bodies. Scattering theory uses the 
refractive index of the particle, the size of the particle and the wave- 
length of the interacting light to calculate the scattering and absorption 
properties of the particle. 

The effects of cloud droplets and aerosol particles on the incoming 
solar radiation can be calculated by using the Mie scattering theory. 

Since aerosol particles occupy a whole range of sizes , the overall 
optical characteristics of an aerosol are obtained by integrating the 
optical properties of each size interval over the whole range of particle 
sizes. The effects of clouds on solar radiation require a complete 
knowledge of the droplet size distribution and concentration ■within the 
cloud (T\7omey, et al. , 1970). Since this data is generally not avail- 
able, we must resort to average cloud conditions. Therefore, we assume 
that the cloud is homogeneous (all particles are the same size) and 
that the concentration of cloud droplets present is directly related to 
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the amount of cloud cover reported. Using these approximations, the 
effects of clouds can he calculated hy using the Mie scattering properties 
of a single sized -water droplet in varying concentrations. Ibis 
technique provides a method of including the effects of aerosols and 
cloud droplets on the atmospheric transmission of sunlight. 

The effects of cloud droplets , aerosol particles and the gases 
present in the atmosphere are combined to produce the overall trans- 
mission properties of the atmosphere. The model developed for this p\jr- 
pose is a modified version of a model developed by Braslau and Dave in 
1972. The model divides the atmosphere into a number of layers. The 
optical properties are considered constant throughout each layer, 
but may change from layer to layer. The model computes the diffuse 
radiation at each level in the atmosphere. The total amount of sun- 
light transmitted throughout the atmosphere to the ground level is 
then calculated by adding the amount of diffuse radiation in the 
downward direction to the attenuated amount of direct radiation. This 
result provides the amount of s-unlight which strikes the earth's 
surface. 

Transmission of atmospheric radiation is calculated using the 
plane-parallel form of the radiative transfer eq.uation. This equation 
calculates the transfer of monocromatic light through a non-homogen eous 
atmosphere and may be -written; 
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- y(T);J(T;y ,^i) (3) 

y = I cos e| where 0' is the angular deviation hebveen the angj e 
of propagation and the local zenith p is defined as upward, 
while -u is downward 
<fi - azimuth angle from the horizon 
. J - is called the soxirce function, it provides the araouiit of . 
radiation added to a specific direction (y, <}>) due to 
■ scattering from all other angles 
T - normal optical depth, it is the svun of optical depth due to 
scattering ancl absorption by molecules and particles 
t(S, M) - optical depth due to Mie scattering (particles) 

t(S, R) - optical depth due to Rayleigh scattering (molecules) 

t(C, M) - optical depth due to Mie absorption (particles) 

W - albedo of single scattering ,• defined as the ratio of optical 

depth due to all types of scattering divided by total 
optical depth 

I - intensity of radiation at optical depth t in the atmosphere,' 
propagating in the direction (p, (j>) 

The source function includes all scattering and absorption proper- 
ties of the mediiim. Essentially it represents the amount of direct 
radiation which is converted to diffuse radiation by interaction with 
the medium. The source function can be represented, .by an equation of 


the form: 
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REPRODUCIBILIT'^ 11; 
ORIGINAL PAGS ?086 


T 

J(t';p,4') = )«F+' 

1 * 

+1 2ti 

Vn* ■// P{T;p,<j>;u' ,<i)’ )' (U) 

-1 0 

I(t;u' ,(j)' )-d.p* 'd(j)’ 

F - Flux at top of atmosphere - • 

t, V,’ ij> ~ have same meaning as "before . 

Vq» <{>q - represent the direction of incoming diffuse radiation 
P the phase function 

The phase function represents the angular distribution of the scattered 
radiation. The normalized scattering function is defined by; 


P(x,v,^;y' ) = T(T)‘M(vi,(i!;p’ ) + 
[1 - T(x)]-R(p,<fr;y',<}.') 


( 5 ) 


T(t) - turbidity factor and is defined as 
- ■ T(t) = At(S,M)/[At(S,M) + At(S,R)] 

M - the normalized Mie scattering phase function 
R - the normalized Rayleigh scattering phase function 
The treatment of the Mie and Rayleigh phase functions is deter- 
mined by the numerical technique .chosen. For use in the iterative 
technique,, the functions are expanded in a Fourier series vhose argu- 
ment is the difference between the azimuth angles of the incident and 
scattered radiation. When this technique is employed it is advantageous 
to express the intensity and source functions in, the form of Fourier ■ 
series also. Since the azimuth integration performed in equation 4 
can be carried out analytically, this reduces the computational 



VI-T 


demand to that of a single integration (Braslau & Dave, 1972). Finally the 
computational 

demand is reduced even further since only the first term of the Fourier' 
series for intensity is needed for the calculation of total raiiation 
flux. Thus, the Mie and,. Rayleigh phase functions may be expressed hs; 



R(Cos8) = ''•P^_^(Cos0) 


M(Cos0) = ^ A. -(Cose) 

k=l ^ 


( 6 )- . 


•- Rayleigh legendre coefficients given by: 

. • ' ■ ‘ 1, Ag = 0, = 1/2 

P - ordinary legendre polynomials 

0 - scattering angle between incident and scattered radiation 
• A^(m).- legendre coefficients for Mie phase function 

The legendre coefficients for the Mie phase function depend only on 
the size distribution concentration and optical properties of the 


aerosol particles. The coefficients are given by: 


^-(m) _ T‘ / Dv(a ,m)-n(r)*dr • (7) ■” 

■ • 4 irB(s,m) r^ ^ 

\ - wavelength of incident radiation 

n(r) “ size distribution of aerosols 
/ 

- unnonnalized legendre coefficients of the series representing" 
the scattering phase fxmetion of a single particle 
. a 7 size parsmeter defined as X = 2tnr/X~ ' 



VI-8 


M - complex refractive index 
B(s,m) - volume scattering coefficient for the medivim 

The volume scattering coefficient can be calculated using the 
efficiency factor for scattering and is given by: 

rg 

B{s,m) = IT' y’ Qs(aj jn) 'n(r) -r^-'dr (8) 

^1 • . ' • 

( 

Q - efficiency factor for scattering 

The unnormalized' legendre coefficients used arise from the expansion 
of 'the con^lex Mie amplitudes by the repeated use of recurrent relationships 
between the derivatives and the products of legendre functions (Deirmendjian 
1969). The method was suggested by Hartel (l970) and first carried out by 
Sekera (l952). The power of this techni<iue lies in the fact that these un- 

normalized coefficients are independent of direction, thus they can be 

} 

computed for the aerosol or cloud particles in advance. We have used a 
modified version of a calculational procedure developed by Dave (19T0) to 
compute these coefficients. This calculation is performed before the 
actual radiation model is operated and are read as input data by the 
radiation model. This technique saved considerable computer time. (See 

Appendix B-1 for program listing.) 

Therefore, the final expressions for the Mie and Rayleigh phase 
functions may be written using the addition theorem of spherical har- 
monies to express 6 in terms of p and as; 
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II 

= I Fj^^^^(p,yM*Cos(n-l)*(.^^-<l.) 


• ( 9 ) 


The only remaining step needed to render the radiative transfer 
eijuation solvable is the addition of two boundary conditions. These 
conditions are defined by assuming that there is no incident diffuse 
radiation striking, the top of the atmosphere and similarly that there is . 
no radiation reflected upward at the base of the atmosphere. These con- 
ditions are expressed as': 

Im(o,-p) = o; Ho' (lO) 


“The effect of ground reflection is included by using a technique- 
developed by Chandrasekhar Cl96o). The intensity calculated by solving the 
previous equations is increased by an amount I^(t, p). The radiation 
reflected by the ground is assumed to be isotropic. The procedure 
■QSgd in-yolves solving the radiative transfer equatzon for the refleci.ed 
ground radiation. The E.T.E. has the form: > . 

y = I«^(T,y) (t,y) 

= h'f P^^^(-r;y,yj^)*e''-- ■" *dyi . 



v/ P^^^(T;y,yn*I»^(x;y^)*dy^ 

-1 • 
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rith boundary conditions given by 


I»°{o;-y) = o 




This -set of equations applies for the case- of Lambertian illumina- 
tion (Ig(-po) = 1). Finally, we have: ^ ’ 

• ±y) = Tg(-iiq)‘I«°(T ; ±p) , , (13) 

^ * X 

The reflectivity R is defined as the ratio of tlie upvard flux given by: 

' 1 2tt 

,Ig(-yo)*y’ f p’dy*d(J) = Tt,*Ig(-po) - (iM, 

o o • ■ 


‘to the total downward flux. The downward flux has contributions, from; 

direct solar radiation TT‘U 0 *F’exp(-T-}j/uo) / (15 

downward flux due to diffuse radiation : • 

1 2ti ■ • • 

J" / l(T-bJ-u,(f>)*dv*d(j) H TTg(T^;-yo) (ic 

o o 

downward flux due to scattered, ground radiation . . 


2 * 7 fIg(-yo)*/l»‘^(T, ;-y)*ydy == TT*Ig(-yo)*S(T^) 


Rearranging the terms in the reflectivity equation ;we. have : 


■ Ig(- o) * 


• R« [yo.F.e-'^b^^° t g(t-b;-yo)3 
• i - R.F(-^ ■)’ 

■ -'b' r 
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The equations presented form, a complete set ■which can he solved to yield 
the flux of solar radiation at any level in the atmosphere* It only remains to 
define the properties of the atmosphere for use in the radiation model. 

Before these equations can actually be used in a .practical model, 
some sin^lifying assumptions should be made to reduce the computational de- 
mand of the model. The effects of including polarization on the total net 
flux calculations has been shown to be insignificant (Adams & Kattawar). 

Thus, polarization has not been included so the scalar form of the R.T.E. can 
be used. Aerosol particles present in the stabilized ground cloud have re- 
fractive indexes which vary considerably with the particle's composition. 

We have assumed an average value of 1,5 for the real part of the refractive 
index. We have also assiomed for this part of our calculation that the 
particles are low absorbing, the imaginary part of the refractive index is 
0.03. Although it is known that aerosol particles in general do absorb 
solar radiation, Braslau and Dave have compared the results of their model 
for the cases of absorbing and non-absorbing aerosols , and their results 
show that only a minor change occurs in the net flux calculations when 
aerosol absorption of less than 0.05 is included (Braslau and Dave, 

1972). 

The radiative transfer equation presented previously is valid for a 
single frequency only. Since most solar energy is emitted between itOOO and 
6000 Angstroms, we have chosen 5000 Angstroms as the wavelength used in 
this model. More accuracy could be obtained, with a drastic increase 
in computational demand, by calculating the groimd flux in each wave- 
length and then integrating the result over the entire solar spectrum. 
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The significant changes in the model required for this calculation are the 
inclusion of the absorption bands for water wapor, atmospheric gases and 
pollutants (See Figure l). Techniques are available to approximate 
these absorption bands (Bergstrom and Viskanta, 19T2). There is no 
reason why the model could not be operated over the entire range, except 
that the co 2 i 5 )Uter time required is excessive. One final simplification 
was utilized, the model is used three times a day at 9:00 a.m, , 12:00 noon 
and 3:00 p.m. local time. The resTilts from these three runs are combined 
to estimate the total radiation received diiring a particular day. These 
simplifications in the model significantly reduce the computational de- 
mand of the model. 

The atmosphere is divided into thirty layers ; the top five layers are 
five kilometers in depth, while the lower twenty-five layers are only ■ 
one kilometer in depth. The stabilized ground cloud is assumed to be 
homogeneoias after T + 3 hours, from the mixing height level down to the 
ground level. The calculation proceeds from the top of the atmosphere 
downward, with the effects of ground reflection added separately. 

A standard atmosphere is used to specify the Rayleigh attenuation 
coefficients and also the Mie attenuation coefficients down to the 
mixing height level. The Rayleightcoeffieients B(s, r) and the Mie 
coefficients B(s, m) are given by McClatchey (1972) for the standard mid- 
latitude atmosphere. These coefficients are used tO' compute the normal 
optical depth and the turbidity for each layer. 


T . 


= B(s,r)i + B(s,m)i 


(19) 



ABSORPTIVITY 
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WAVELENGTH (microns) 


FIGURE 1 

Absorption spectra for H^O, CO^, 0^, K 2 O I'ltid the atmosphere. • 
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(t) = AT^(s,in)/CAT^(s,m) + AT^(s,r)H (19) 

The legendre coefficients for each layer are calculated using 
equation 7, where the size distribution is assumed to be constant after 
T + 3 hours. Thus, the optical properties of all thirty layers are com- 
pletely determined for the standard atmosphere. A value of .3 was used 
for the albedo. This value is read into the model as data so that it 
can be changed quickly if desired. 

If clouds are present during the calculation period, their effects 
are included at the three kilometer level and below. The cloud drop- 
lets are all assumed to be average sized stratus cloud droplets 
(h microns in radius ) (Mason 19Tl) • The concentration of droplets is determinec 
by converting local climatological cloud cover data to tenths of sky 
cover. The cloud ceiling reported by local weather stations is used 
to determine the thickness of the cloud layer. Since the scattering 
efficiency factor for 4 micron sized water droplets was previously 
calculated (Appendix B-1) , the attenuation coefficient for cloud 
layers B(s,c) can be computed using equation 8. Therefore, the total 
attenuation coefficient for the layers is just the sum of B(s,m) and 
B(s,c). Finally, the legendre coefficients for the layer are simply 
the standard coefficients plus the coefficients computed for the cloud 
layer using equation 7. This determines all the optical properties of 
the cloud layer needed for the flux calculation. 
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The mixing depth determines the height in the atmosphere pene- 
trated by the stabilized ground cloud. The aerosol mass concentration 
of the stabilized ground cloud is determined by the assumed dispersion. 
The mass concentration is converted to a number concentration of 
particles between the reported size interval. The total number of ■ 
particles in each .1 micron size interval is computed. The results of 
this computation are used in equation 8 to determine the new attenuation 
coefficient for the layers. This procedure is similar to the procedure 
used to determine the attenuation coefficient for layers of cloud drop- 
lets. Finally, the legendre coefficients are computed using equation 7. 
The new values for the legendre coefficients and the attenuation co- 
efficients replace the standard coefficients for all levels below the 
mixing height. 

Recalling equation 4, it is seen that the first term of the 
source function is the contribution due to scattering of direct solar 
radiation to the source function. This is the primary scattering term. 
The iterative technique approximates each successive order of scatter- 
ing by using the previous result to obtain each higher order scattering 
term. Thus, the primary scattering is approximated by: 

= Js*F*exp(-ri/po) *P2^(T^;^y ,-po) (20) 

where P is given by equation 5. 
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The successive approximation of the intensity function is then 
given hy: 

^ ^ (Tj,_^;-v) *exp{.-AT^/v) .+ 

' ’ *[l “ exp(-ATj^/y)l ■ 

for the upvard intensity and similarly for the do-wnvard intensity. 
Successive orders of scattering are included in the so\irce function by 
approximating equation 4 as: 

. ' ;±y) + h’ 

X J. * 

J'Pl(T^;±y 3 -uM*Ii^^^^(Ti;-n^)*du^ + (22) 

o 



o 


Ja is calc\ilated using a sub-program 'called SOURCE for listing see 
Appendix B-3.- This approximation is substituted in equation 21 to 

t 

obtain the intensity function vhich includes the next order of scatter- 
ing. This new intensity is then inserted into equation 22 to calculate 
the new source function which now includes scattering of the next order 
The procedure is continued until the desired agreement between Im and 
im+1 is achieved. The criteria used for this model terminates the 
procedure when; 
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■ (lm+1 - Im) ^ .001»Im (23) 

. The effects of ground reflection are included using a procedui'e en- 
tirely analagous to that presented above.’’ The only significant differ- 
ence being that calls to SOURCE are replaced, vith calls to sub-prograjn ' ' 
GROUND. (For listing see Appendix B-4. ) The procedure is applied 
only once for the gro\md reflected radiation since it_ is assximed to be 
reflected isotropically at the ground. 

■ The model finally integrates the intensity function, which now in- . 
eludes all effects due to ground level reflection, through- all do’v.Tiward 
angles to obtain the net downward flux. This is accomplished using the 
■flux equation given by; ■ 

F^(t) = 2’TT'y I^^^(T;-p)*y‘dp , (24) 

•o' 

. This flux value represents the total amount of diffuse radiation ' 
reaching the ground level. The directly transmitted radiation is com- 
puted and added to the diffused using: 

■ . Flux = Flux + lQ’exp{-T-jj/yo) . • (2?) ' 

s 

As mentioned previously the model is used three times- a day to- 
obtain the fin^ value which gives the total flux for the entire, day. 

A listing of the main program is provided in Appendix B-2. 
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The model was tested at the city of Albany since measured values 
of solar radiation are also available. The results of the week-long 
model validation run are presented in Figure VI-2, A. second plot of the 
measured values of solar insolation for the same period is also included 
in the graph to facilitate the comparison. In Figure VI-3, we have plotted 
the amount of clouds present and also the cloud celling for the same time 

. V 

period. From a study of the graphs, it is seen that the calculated values , 
and the measured values of solar insolation agree quite closely for the 
period . 

The model was then used to demonstrate the effect of the stabilized 
ground cloud on solar energy transmission for the Florida latitude. Cloud 
cover was set equal to zero for these calculations. 

Figure 4 shows the overall results of percent decrease of ground 
level solar radiation versus aerosol mass concentration within the 
stabilized ground cloud. We have carried out the complete radiation 
transfer calculation in an atmosphere clear of water clouds. The limit 
of vertical mixing was assumed to be 850 meters and 1550 m. Below the 
mixing height, the stabilized ground cloud is homogeneously mixed in the 
vertical direction. Figure 4 constitutes the basis for our impact assess- 
ment with respect to the radiation flux received at ground level. 
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RESULTS 

1. Any effects resulting from the absorption of solar radiation by 
the stabilized ground cloud (deposition of heat energy in the atmosphere) 
are negligible, if the absorption coefficient of the particles within 
the S.G.C. is less than 0.03. 

2. The scattering of solar radiation by the stabilized ground 
cloud resvilting in a decrease of solar energy received at the ground 
amounts to maximal 15 ^. 

Under those high scattering conditions, the S.G.C. •wo\fLd have to be 
dispersed homogeneously between ground and 1550 m with an average mass 
concentration of 15 O pg/m^ which could occur at times up to T + 3 hrs. 
Eor times greater T + 3 hrs. , we anticipate the decrease of solar radia- 
tion at grotind to be less than 15^ under all meteorological conditions . 
This change of heat energy deposition at the ground is not sufficient 
to cause significant changes in the vertical mixing of air causing such 
effects as cumulus cloud formation, etc. We, therefore, conclude that 
there will be no detectable impact on ’’weather" resulting from the 
interaction of the stabilized ground cloud with solar radiation. 

However, we anticipate a reduction in "visibility". Since there will 
soon be a secondary air quality standard with respect to "visibility" , 
some attention should be devoted to the radiation-interaction problem. 
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ees'eaech recommemdatioms 

For more precise calculations , the real and imaginary part of the 
extinction coefficient of the S.G.C. aerosol must be knoTOi. If the 
aerosol absorption coefficient exceeds 0 . 05 j then the problem must 
be reassessed. 
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FIGURE VI-1 


Comparison of measured and calculated values of solar insolation for 
six“day period. 




FIGURE 'VI-2 ■ 


Cloud cover and cloud ceiling height. for the six-day modeling period. 
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'iraNOKlIfiLIZED' LEGENURE XCFEFFXCIEKTb" PKaGESIT' LISTING- 


INTEGER PtQ 

REAL TA(4) ♦T8(2) *TC{E) *T(4) »COEF( 100) 

. real NOC,N!OB»KOC,KOB ... 

COMPLEX BCAp{ 500) «AD0D<500) 

COMPLEX RF »RRF ,RRFX»WM1 »FnA« FNB »'TC 1 oTC2» WFN(3) sACAPiSOO) 
COMMO.NP » G # X « AOOD . ACAP » CKR 

EQUIVALENCE <WFN( 1 ) ,TA( 1 ) ) » (FNA»TB( 1 > ) . (FNBtTCC X > J 
. REAO(5*5)RFR,RFl sX 
S FORMAT t3FlO«2) 

DO 300 LIN=1«6 
X=* I^LIN 

2 RFR=1.5 

3 X=2.*3« 1416*-X/*5145 
RF=CMPLX(RFRt-RFl ) 

RX=1»/X 
RRF=1,/RF • 

RRFX=RRF^t-RX 

T( I ) =2^tX+10« 

NMXl = i • iO^TI 1 ) . 

NMX2=T<1) 

IFlNMXlcGTolSO) GO TO 22 

NMX 1=150 

NMX2=135 

d.Z ACAP<NMX1 + 1 ) = (0» *0* >. 

THE FIRST LOOP COMPUTES THE MIE FUNCTIONA AND B USING A 
DOWNWARD RECURSION TECHNIQUE 
DO 23 N= I tNMXi 
NN=NMX1-N+1 

ACAP(NN) = <NN+1 )^4^RRFX-1 a0/( (NN+X )*RRFX+ACAP{NN+1 ) ) 

23 CONTINUE 

T ( 1 ) =COS.( X > 

T(2)=SIN(X) 

WM1=CMPLX(T{ 1 ) ♦'*T(2> ) 

V\;FN ( 1 ) =CMPLX ( T (-2 ) ^ T ( 1) > 

WFN { 2 ) =RX*WFN { 1 ) -WM I 
TCI=ACAP{ I )*RRF+RX 
TC2=ACAP( 1 )*RF+RX . 

ACAPU ) = (TCI*TA(3)-TA{ 1 ) ) /< TC 1 -»WFN ( 2) -WFN{ I > ) 

8CAP( i ) =CTC2*TA(3)-TA( 1 ) ) /< TC2*WFN ( 2) “WFN { 1 ) ) 

FNA=ACAP( X ) 

FNBsBCAPU) 

. TC1)=3.0 

QEXT = (TaC 1 >.+ TC{ I > )*T< 1 ) 

GSCAT=<TB< 1 )-»-5t2+TB«2)**2+TC{ 1 ) *5-2+TC(2)**2)*T( X) 

N«2 

65 TCl>=2«N~l 
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T(2>=N-1 
T(3)=E*N+1 
WMl=WFNtl) 

WFN( 1)=WFN(2) 

WFN(2) =T{ t )*RX*WFN( 1 )~WMl 
TC 1 = AC AP ( N } *RRF+N*RX 
TC2=ACAP ( N ) »RF+N*RX 

ACAP (N)= (TC1*TA<3)-TA< 1 ) > / < TC 1*WFN < 2 > ~wFN ( 1 ) ) 
BCAP(N)={TC2*TA<3)-TAU ) ) / CTC2»WFN ( 2 ) -WFN U ) ) 

FNA=ACAP(N> 

FnB=BCAP(N) 

T( 4)=TB( 1 >**2+TB( 2)^*2+TC{ 1 J-»*2+TC(2 3*-»-2 
QEXT =QEXT+ t T ( 3 ) *TB ( 1 > +TC < 1 ) ) 

QsCAT=QSCAT+T{3')*T (4> 

N = N+l 

IF<N«LE.NMX2)G0' TO 65. 

WRITE(6.50) T(4) 

50 FORMAT (@ TO BAO@*E20o6> 

100 T( 1 )=2.0*RX**2 
NUMT=N 

QEXT=QEXT*T( 1 ) 

QSCAT=QSCAT*T{ 1 ) 

WR 1 TE C 6 , 6 ) QEXT « QSCAT 
6 F0RMAT{2E20*8> 

N=N/2*2 
N=N+2 
K = N+1 

C- THE ODD AND EVEN SUMS OF A AND B ARE COMPUTED FOR LATER USl 

AD0D(N)=0*0 
AD0D<K)=0.0 
110 N=N~2 
K=N-1 

NOC=l./N+< 1 o/lN+l ) > 

K0C=1 ./K+( 1 ,/(K+l ) ) 

KOB=l ,/(K+l > + ( 1 o/(K+2) ) 

N0B=1,/(N+1 )+{ 1 */(N+2) > 

ADOO(N)=BCAP(N)*NOC-NOB*ACAP(N+1 )+AD00(N+2) 

ADOO < K ) =BC AP ( K ) -»KOC-KOB*ACAP < K+ 1 } +AOOO ( K+2 5 
IF(N.N£.2) GO TO 110 
DO 104 K=1 1 12 
I=K*10 

WRITE<6,102) ACAPU ) ♦BCAPC I ) , ADODC'l ) 

102 format ( 6 E 2 O 98 ) 

104 CONTINUE 
K=0 

THE MAIN program CALCULATES AN SERIES' APPROXIMATION 
TO the UNNORMALIZED LEGRENDRE COEFFICIENTS NEEDED IN ■ 

the rte model 
ktot=numt -2 
TMO=oO 
BKT?0*0 

105 K=K+I 
TMO=OoO 
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8KT=0.0 
KC= <K+1 )/2*2 
IF(KC.EQ.K>“ GO TO 175 
, KPRIMs(}<-l >/2 
M-KPRIM ■ 

iF<K*GT*l > GO TO 1 IS 
A10=2. ■ 

. All=AlO- 
BI0=1. 

P = M+l 
Q«p 

CALL. LSUM 
BKT=810*CKR 
GO TO T60 
IIS 1=0 

A10=4*<K-1 )*(K-2>*A10/< (2-»K-1 )^{2*K-3> ) 

A11=A10 

BlO=(i<-2)**2*BlO/CK-l )**2 
811=810 
Q=M+I+1 
P=M“I+1 
CALL LSUM 
BKT=B1 1*CKR+BKT 
IF ( I •EQ-KPRIM) GO TO160 
GO TO 150 
lAO M=M+1 
1=0 

A1 1 = {2^M-K)*(2-5^M-1+K)*A1 1/1 1 2*M+K)*(2*M-K+1 ), ) 

P=M+1 

Q=P 

CALL LSUM 
BKT=61 O^CKR+BKT 
811=810 

IF( I •EQ.KPRIM) GO TO160 
ISO 1=1+1 

Bi 1=(K-2*I + 1 )*(K+2*I-*2)*B11/1 (K-2*I )-»(K+2*l~l ) ) 
P=M-I+1 
Q=M+1+1 
CALL LSUM 
BKT=BU*2*CKR+BKT 
IF ( I .LT.KPRIM) GO TOISO 
160 TMO=Al 1*BKT+TM0 

IF(M.LT»KT0T) GO T0140 
COEF(K) = U<-*5)*TM0' 

COEF ( K ) = A8S < COEF IK)) 

IFIK.LT.NUMT) GO • TO 105 
GO TO 220 
175 KPRIM=(K-2)/2 
M=KPRIM 

IF(K.GTp2)GO to 185 
AD 1 0=.4 « /3 • 

AOl l=A010 
8010=1*0/2»0 
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185 . 


ISO 


195 


200 


220 

230 

240 

250 

300. 


BD 1 1 =BO 1 O - 
P=M+1 
Q=P+l . 

CALL. LSUM 
BK.T = BOlO*Ci<R*2 
GO TO 200 
1=0 

AO10=4*(K-l >*(K-2)*ADlO/t (2*K-1 >*(2-»K-3) ) 

AD 1 1 =AD 1 0 

BOlO=<K-l )*(K~3>*BD10/<K*<K-2) ) 

BD11=8010 . ■ ■ 

P = M + i 
Q=P+1 

Bi<T=BO 1 1 «-CKR-3t2+BKT 
CALL LSUM 

IF tl oEQ.KPRIM) GO TO 200 
GO TO 195 
M=M+1 
1=0 

A01l=(2*M-)<+l)-5t{2*M-!-K)*ADll/'( (2*M+K+1 )*<2*M-K+2) > 

P=M~1+1 

Q=M+I+2 

CALL LSUM 

BKT=0D 1 0*CKR*2+BKt 
BDI 1=8D10 

IF ( I .EQ*KPRIM) GO TO 200 
1=1 + 1 

BDI 1=( 2+I+K-l )*(2+I-K)-»BDl 1/ ( (2*I-K+1 >*(2*I+K)> 

P=M-I+1 

Q=M+I+2 

CALL LSUM 

BKT=8D 1 1 +2^^CKR+8KT 
IF ( 1 .LT.KPRIM) GO T0195 
.TM0=A01 1*SKT+TmO 
IF (M«LT«KTOT > GO T0180 
COEF(K)= (K-*5)*TMo 
COEF<K)=ABS(CO eF(K) ) 

IF (K*LT*NUMT ) GO TO 105 
WRITE{6»230) (COEF(K) ,K=1 »NUMT) 

FORMAT <5E 14,4) 

WRITE(7,250> (COEF(K) ,K=1 ,NUMT) 

FORMATtSElO.A) . 

UINA=JINA+1 
END • 


SUBROUTINE LSUM 

LSUM COMPUTES THE INNER MOST SERIES OF THE APPROXIMATION 
FO USE IN THE MAIN PROGRAM 

integer P»Q 



COMPLEX ADC©(500) »ACAP(500) 

COMPLEX AiAl^D»E 
COMMCfxJP , Q , K , ADOD » ACAP « CKR 
IF (P.EQ. 1 ) GO TO 10 
A=ACAP(P-1 ) 

Al=ACAP(Q“l } 

5 D=(2*P-1 )*<P~1 )*A/P+(2*P-1 }*AOOO<P> 
E=(2*Q-1 )-5-(Q-l )*Ai/Q+<E*Q-1 )^AD0D(Q> 
E=CONJGIE) 

CKR=REAL(D*E) 

RETURN 

10 IFIP.EQ.Q) AI=CMPLX(0»0*Qo 01 
A=CMPLX(O.O»0e0) 

GO TO 5 
END 

@END ■ • , 

@ELT*0IL D*0- 

1.33341 0.0 4.5 

©END 

@FIN 


000232 
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APPENDIX B-2 

R^IATIVE TRMSFER MODEL — MAIN PROG^ LISTING 
COMPZL-ER (DIAG=3> ” 

real int(3o«20) * Into 

DIMENSION CL0U(3) » CLOUCI 3 5 a CCOEF ( 34 ) *eETA UO ) » BETR ( 10) • * "■■■■■ 
X BTA< 10) »QSCAT<7) ,C0EF(6<i34) »ACO( 1 I » 34 ) »VAL< lO) « 

X TURB( 10) .TAU(31 ) »C0NCC3) ^AER( lO) »FAC(6) . 

X PHASE <30. 20* 1 1 ) »GNT<30.20) .SOUR <30, 20) *J)SqUR ( 30 ♦ 20 ) *HMXT{3> ’ 
X.FA<6) ,BTR< 10 j .CA<6) 

,COMMONPHASE, INT.GNT.LI * J I iDELT. I tOSOUR.SOUR 
•DATA /FA/. 1984, .aElS. . 1462, .1091 , «0932» •072S/'- 
READ <5, 26) BETA 
READ(5,26)BETR 
26 FORMAT (5F 15, 10) 

READ <5,39 )AER 
39 FORMAT (5E1 0.4) 

REA0C5.2B)QSCAT 
28 F0RMAT(7Fr0o8> 

• !rEAD<5.22)CC0EF 
2.2 format (8E 10.4) 

00 30 1=1,6 

30 REA0(5.22) <CO£F< I ,K> ,K=1 ,34) 

BET0=1,25E“6 
bet 1=1 ,64E-S 
TAU (31 ) =0, 

AMDA=.5i45 

FLUX=295, 

- CR=4, 

NUM=34 

ELEV=.247 

CORR=lo-ELEV 

PI =3 ,14 16 

NUMT=0 

CBET=0, 

DO 46 0=1,20 
I'F< J.LE, 10)0NT< 1 , J)=0* 

46 IF ( J.GT. 10)ONT<30»'J^ ♦ 

READ(5, 17) AUBED 

17 FORM AT (Flo, 2) ^'EPEODUCIBILITY -OF THE 

0AS=160, '^IGINAL PAGE IS POO^' 

c=o. 

DR=,1 

co=o 

. DO 47 1=1, 100 
R=.1*I 
C=0R/R«-*4-fC 
IF ( I ,ET * 50 )CO=DR/R+CO 
IF< I eLE^eiFACt 1 )=0R/R*-»4 
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47 IF( I .l_E®6)FA( I )=DR/R 
C=1 •/€ 

CO=i ./CO 

REAO(5.53)Cl *Ca*C3 

50 DAS=DAS+1 
JIN = 0 

53 FORMAT(3F10.2) 

JIN0=0 

51 IF( JINO.EQ. 1 )O0 TO 55 

REAO(5.20)CLOU 
DO 52 5=1 » 3 

IF(CLOU( I ) .LT. .aiCLOUC i )=CLOU( I )*C1 
IF(CLOU( I > .GT* .2)CL0Ua )=CLOU< I )*C2 

52 lF(CLOU< 1 ) *GT..8)CL0UC I )=CUOU( I >*C3 
20 FORMAT (3F1 0*2) 

READ(B.20)CL0UC 
READ(5«20)HMXT 
REAO(5*20>CONC 
55 OlN=JlN+l 
NUMT=NUMT+1 

CU0UC( J I N ) aCLOUC (JIN) /3280.8+ELEV 
. I CE I L=a_0UC ( J I N ) + * 5 
• IHMXT=HMXT( JIN)+.75 
IF (CLOU< JIN) .LT. -01 > GO TO 200 
CBET=2.E-1 1*CL0U( JIN) 

CSET=P I *OSCAT { 7 ) *CR**2*C8ET 
C ATT^NUATldN COEFFICIENTS DuE TO THE PRESENCE OF CLOUOS 

C ■ ARE COMPUTED USlNG THE SCATTERIN PROPERTIES OF 4 MICRON DR 
op 130 K=1,NUM 

ACO< 1 1 ♦K)=AMDA.-»*H*CCOEF(K)*2.E-10*CLOU( jiN)/(PI*CBET) ’ 

180 -CONTINUE • ' 

CSET=CBET*1 .E9 

C SCATTERIN COEFFICIENTS DUE TO AEROSOLS ARE COMPUTED IN 

C THIS LOOP. 

200 DO 250 L=1 ♦ 10 , 

' IF(CLOU( JIN).LT.«01 >CBET=0.0 
BET=BETA(L) 

BER=BETR (L) 

IF(L.GT.4>BET=BET-BETACL-1)- 
IF(L.GT.4)6ER=3ER-BETR(L-1 ) 

BTA(L)=BET 

BTR(L)=SER 

■ i F( ( lO-IHMXT ) .GE*L) GO TO 235 
DO 215 I=l«6 
R=.1*I 

CON=CONC< JIN)*PA( I )*C0*3.E-12/(4.*P1*R**3) • 

IF{ { lO-lCEIL) .LT*L.A nO. ICEIL*N£.0) GO TO 214 
IF(CLOU( JIN) .GT«0» )C0N=C0N-» 1*C0N*CL0U( JIN) 

214 VAL( I )=0SCAT( I )*C0N*R*-if2*PI 

215 CA( I )=CON 

8ET=.5/18.*(vAL( 1 >+4-»VAL(2)+2*VAL<3)+4*vAL(4) 

X +2-»VAL(5)+VAL{6> ) 

. DO 230 K=1«NUM 
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DO 220 1 = 1 »6 
R = • 1 [ 

220 VAL(I)=COEF{I.K>*CA<I> 

AGO (LtK)=AMDA**2*.5/ (PI*BET*18, )* 

X < VAL( 1 )+-4*VAL(2>+2*VAL(3)+4*VAL(4)+2*VAL<5)+VAL(6) ) 

IF (( 10- 1 CEIL) .GE.G) AC0<L«K)=AC0{L»K )+ACO( 11 »K) 

230 CONTINUE 
GO TO 246 
235 8ET=8ET*1 *E-9 
DO 24S K=1 »NUM 
DO 240 1=1 «6 . . 

R=*l«-I 

• CON = AER ( L ) *FAC ( I ) *C 
IF ( ( 10-ICElL) •GE»L.0R.L.LE.6 >GO TO 240 
IF(CLOU( JIN) .GT^O* )CON=CON-o 1 -*C0N5<-CL0U ( J IN ) 

240. VAUt I >=COEF( I ♦K)*CON 

AGO ( L, K ) = AMDA**2* . 5/ ( P I *8ET*- 1 a . ) 

X*(VAL( 1 )-t-4*VAL(2>+2*VAL(3)+4*VAL(4)+2*VAL{5)+VAE(6) ) 

' IF( ( 10-ICElL) oGE.L. 0R.L.LE*6) GO TO 245 
AGO(L,K) iACO(L,K)+ACOU 1 »K) 

245 GONTINUE 
■ 246 GONTINUE 

BET=BET*1 *E9 
BTA<L)=8ET 

TURS ( L) = ( BET+GBET ) / < BET+CBET+BER ) 

IF I ( 10-IGEIL) *LT.L.0R«L.LE,6 ) TURB ( L ) =BET/ ( BET+8ER > 

250 GONTINUE 

40 FORMAT (5E 10.4) 

THE AZIMUTHAL ANGLE OF THE SN IS- COMPUTED USING THE- NUMBER 
OF DAYS SINCE THE SUMMER SOLSTICE 
RND=2,*P 1/360. 

AT=43.*RND 

DEC=DAS*360./36S .2563^RnO 
DEC = 23 • 4S*RND*C0S < DEC ) 

H=ABSt45.*RND*(NUMT-2) ) 

PMUO= ( S I N ( AT ) *S I N ( DEC ) +COS ( AT ) *COS < DEC ) *COS ( H ) ) 

PMU0=-1 o*ABS(PMUOl 
WR ITE(6»41 )PMUO 

41 FORMAT <E1 0.4 ANGLE®) 

IFIDAS.LT. 162. .and. JIN. EQ. 1 )HUX=FLUX 
FLUX=HUX*.225 

IF < U IN.EQ. 2)FLUX=HUX*.55 
WRITEI6. 191 )FLUX 
19l' FORMAT!® FLUX = '®.F4.0) 

■ T = 0. 

the LOOP COMPuTESTHE SCATTER in PHASE FUNCTION FOR USE 
IN LATER CALCULATIONS 
DO 350 L=1 .30 
IF(L.GT.24) GO TO 270 
IFIL.GT.5) -GO TO 256 
N=1 

8ET=L*8TA ( N ) +BET0 
BER=L*BTR ( N ) +8ET 1 



T AU { L) = { BEt+BEP ) +T 
T=TAU(L) 

8ET2=BET 
BET3=8ER 
GO TO 365 

256 IFtL.GT.U) GO TO 250 
N = 2 

8ET=8TA(N>+BET2 
8ER=8TRCN>+BET3 
BET2=BET 
BET3=BER 
GO TO 265 

258' IF(L.GT. 16) GO TO 260 
N = 3 

BET=BTA<N)+BET2 
BER=BTR{N)+BET3 
8ET3=BER 
BET2=BET 
GO TO 265 
260 N=4 

BET=BTA(N)+BET2 

BER=BTR(N>+BET3 

BET2=8ET 

8ET3=BER 

265 TAU(L)=8ET+BER+T ‘ 

T=TAU(L) 

GO TO 300 
270 N=L~20 

TAU { L) -BTA ( N ) +8TR 1 N ) +C8ET+T 
I F t N .LE • 6 ) TAU 1 L ) =T AU ( L ) -C8ET 

IF (L.EQ»30 )TAU'( L> =bTA(N) *C0RR+BTRC N) *CORR+CBET*CORR+T 

IF ( <30-L5 ,GE. ICEIL) GO TO 290 
I F ( L .EQ . 30 ) CBET=CBET*CORR 
T AU ( L ) =T AU { L ) -CBET 
290 T=TAU(L) 

300 continue 

DO 350 J=1 « ii 
PMUP=*-J*. 1 

I F ( U ,E0 .11) PMUP=-PMUO 
278 DO 350 1 “I .20 
PMU=-I*.l ■ 

1F(I.GT«10> PMU=( I-10)*« I 
PMU2=l. . , 

PPMU2= 1 . 

FONEM= ACO ( N , 1) *PMU2*PPMu2 

PPMU3=PMUP 

PMU3=PMU 

DO 320 K=2.NUM 

FONEM= ACO ( N . K ) *PMU3*PPMU3+P0NEM 
PPMU1=PPMU2 
PMU1=PMU2 
PPMU2=PPMU3 
• PMU2=PMU3 
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PMU3=2.*PMU^fPMu2-PMUl-(PMU*PMU2-PMUi j/CK+l ) 

PPMU3= 2 • *PMljP*PP'''tU2-PPMU 1 - < PMUP*PPMU2-PPMU V ) / < Ki- j ) 
320' CONT INUE 

FONER=3 • * < 3 • *-PMU**2-PMUP**2+3*PMU**2*PMUP*«2 ) /8 • , 
PHASEtUf I ♦ J)=TURB(N)*F0.NEM+C 1 <,~TURB(N) ) *FONER 

350 PHASE<L. I ♦ J)=ABS(PHASE(L» I »J) ) 

CP = 0 

CPO=0* 

DO 352 l_=l ♦30- 
DO 351 J=1 t20 
DO 351 N=1 ♦ 10 

351 CP=CP+PHASE{L* J«N) 

CP=1*/CP 

DO 352 J=1 ^20 
DO 352 N= 1 * 1 0 

352 PHASECL, J»N)=PHASE(L»J^N)*CP 
DO 354 L-1^30- 

DO 353 J=»1^20 

353 CPO=CPO+PHASE(L« Jo 1 n 
CP0=1 ./CPO 

DO 354 J=lo20 

354 PHASE ( L ♦ 3 ♦ 1 i ) -PHASE ( L » J ♦ 1 I ) *CPO 
45 F0RMAT{5E10«4,@TATATATA@) 

WRIT£(6«45) (TAU(I-) ♦L=25,30) 

THE FIRST approximation TO THE INTENSITY FUNCTION IS 
CALCULATED INCLUDING ONLY FIRST ORDER SCATTERING 
DO 380 L=lf30 
DO 380 J=lolO 
DO 380 1=1 »2 
0TA=TAU(L) 

LB=-1 

JI=J 

LI=L 

D=U 

PMU=-* 1*0 

V=-l .*ABS(DTA/PMU) 

V=EXP(V) 

IF ( I ,EQ. 1 >G0 TO 356 

LI=31-L 

LS= I 

J I = J+1 0 

PMU=-PMU 

DTAs-l .*( A8S(TAU(LI ) -TAU ( L I -f LB ) ) ) 

V = -l .*ABS(OTA/PMU> ‘ ' 

V=EXP< V) 

356 IF(L.NE.1)G0 TO 370 
InTCLI 0 JI )=0o 
GO TO 380 

370 DTA = -1 .*ABS(TAU(>-I )-TAUCL1+L8> ) 

V=-l •*ABS(DTA/PMU) 

V=EXP(V) 

OELT=TAU<30)~TAU1L1 ) 

W=-i*A8S(TAU(LI+l >/PMUO) 
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OSOUR ( L I » J i ) =FLUX*eXP ( W ) 

X *PHASe<LI «JI» 11 > 

INTCL! , JI )=INT(LI+LB» JI )*V+OSOUR(LI * Jl )* ( 1 .-V) 

CALL GROCNO 

GNT(LI tJI >=GNT<LI +LB» JI )*V+SOUR(LI ♦ JI>*( I .r-V) 

380 CONTINUE • . . 

400 DO SDO J=l,10 

SUCCESSIVE ORDERS OF SCATTERING ARE INCLUDED IN THE 

intensity calculation until comparison criteria is met 

DO 500 I = 1 *S 
DO 430 L=1 »30 • • 

DTA=TAU<L> 

L8=-l ’ ■ . 

JI=J 
LI=L ■ 

PMU=-J*9l 

IF( I .EQ* 1 >G0 TO 450 
LI=31tL 
•L8=l 
JI=J+10 
PMU=-PMU 

450 IFI.LpNEoI) GO TO 475 
GO TO 480 
475 call source ' 

DTA = -1 «*( ABSITAUILI )-TAU(Ll+LB) > )' 

V=-l .*ABS(DTA/PMU) 

V=EXP(V> 

IF! JI,EQ.lO) INT0=INT(LI » Jl > 

1NT(LI « JI ) =SOUR (LI , JI ( U-v )+ 1 NT (L I +LB , J I.) 

480 CONTItslUE 
500 CONTINUE 

WR 1 TE { 6 ♦ 482 ) I NTO ♦ I NT ( 30 , 1 0 > 

482 FORMAT(2E10»4«©INTOINTO©> 

IF{ A8s( InT(30i 10) -INTO) .LE« » 001 «-INT, 0 )GO TO 510 
• GO TO 400 ' . 

THE EFFECTS OF GROUND REFLECT I ON ■ ARE INCLUDED IN THE FJ NAL 
CALCULATION OF THE DIFFUSB COMPONENT OF THE GROUND LEVtL 1 
510 DO 550 J=1 « 10 

PMU=-,1*J • 

VAL( I )=INT(30* J)*-PMU 
550 S0UR(2* J)=GNT(30, J)*-PMU 
■ I = l 

540 TEGRAL^^g/SO^-s-tVALt 1 > +4* VAL ( 2 ) +2 *VAl ( 3 ) +4»VAL ( 4 ) +2*VAL ( 5 ) 

X +4*VAL(6)+2*VAL(7)+4*VAL{8>+2-*VAL{9)+VAL( 10) ) 

555 SOURd « 1 )=2.*TEGRAL 
1 = 1 + 1 

IF ( I •EQ»3>G0 TO 561 
DO 560 J=1.10 

560 VAL{ J)=S0UR(2» J> 

561 CONTINUE 

IF{ I «EQ»2) GO TO 540 

REFL= ALBED*(PMU0*FLUX+EXP( TAUIGOl/'PMUO+SOURta* U } 

X /U-AL8ED*S0UR(l»i>) 
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DO 580 J=1 i 10 
6NT ( 30 » O ) =REFL*GNT { 30 t O > 

INT (30» J)=INT<30» J )+GNT( 30» J ) 

580 VAL< J)=INT(30» J)** 1*J 

TEGRAL= . 9/30 * *■ ( VAL ( 1 ) +4*VAL ( 2 ) +2*VAL ( 3 ) +4*VAL ( 4 ) +2«-VAL ( 5 
X +4*VAL(6)+2^tVAL(7)-h4*VAL{8)+2*VAL(9>+VAU(l0) > 

INT(30» 1 )=TEGRAL*2.*PI 
FLU=-1 *ABS ( T AU ( 3o > /PMUO ) 

INT(30«1 )=lMT<30i 1 ) FFLUX»EXP (FLU ) 

WRITE{6«a0) INT(30»1) 

80 FORMAK 1E14«5) . , 

TIN sslNT(30i 1 > + TlN 
lF(JINeNE«3) GO TO 55 
WR1TE{6«585>TIN 

585 FORMAT <@ SOLAR FLux©»3XiF6o2 ) 

TIN =0 . 

numt=o 

JIN=0 

IF(DAS*LT. 166)60 TO 50 
END 


' 000338 
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RADIATIVE TRANSFER MODEL — SUB-PROGRAM SOURCE LISTING 


SUBROUTINE SOURCE 

SOURCE COMPUTES SUCCESSIVE APPROXIMATIONS TO THE SOURC . 
FUNCTIONWHICH include higher ORDER SCATTERING 
real INT(30»20) 

DIMENSION VALC 10) »0S0UR(30»20) 

X » phase ( 30 . 20 » 1 1 ) » GNT ( 30 » 20 ) ♦ SOUR <30*20) 

COMMONPHASE* I NT*GnT * L *K » CELT *' I I *OSOUR*SOUR 
PI=3*1416 
■ S0UR<L*K)=0 
J0=i 
JI=JO+9 
KI=K 

DO IS 1=1*2 
5 DO 10 J=UO*JI 

10 VAL( J-JO+1 )=PHASE(L»KI )^INT<Li J) 

TEGRAL= • 45/30* * ( VAL ( 1 > +4*VAL ( 2 ) +2*VAL < 3 > +4*VAL < 4 ) 

X +2*VAL(5)+4-»VAL<6)+2*VAL(7)+4-»VAL<8)+2-i^-VALC9)+VAL< 10) ) 
JO=ll 
JI=JO+9 

IF( I UEQ.2)KI=K-10 
IF( II*EQ« 1 )KI=K+10 

1 5 SOUR ( L * K ) =OSOUR ( L , K ) + TEGR AL+SOUR < L « K > 

RETURN 
. END 


■ 000035 
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APPENDIX B-4 ■ . • 

RADIATIVE TRANSFER MODEL ' — SUB-PROGRAM GROUND LISTING 

» 

SUBROUTINE GROUND 

SUBROUTINE COMPUTES THE SOURCE FUNCTION NEEDED TO 
■ INCLUDE THE EFFECTS OF GROUND LEVEL REFLECTION IN THE 
INTENSITY CALCULATION 
real I NT (30*20) 

DIMENSION VAL( 20) *OSOUR(30*20) 

X* PHASE (30 *20* 1 I ) *GNTC30*2O) « SOUR (30 *20} 

■ COMMONPHASE* INT»GNT*L»Ki DELT 1 I I »OSOUR»SOUR 
LI=L • • 

LB=-1 

IF( I I.EQ.2)LB=1 
S0UR(L*K)=0 
TEGRAL=0 

DO 1 0 J=1 *20 

JI = J 

KI=< ■ ■ ■ - • . 

■ IF( J.LEo 10) GO TO 8 
IF(K.LT. 10>KI=K+10 
IF(K.GT. 10)KI=K-10 
JI=J-10 

8 CONT INUE ■ ....... 

VAL{ J) =PHASE(L,KI »U1 ) *GNT(L1 +LB* J ) 

10 CONTINUE 

TEGRAL= . 45/60 . * < VAL ( 1 ) +4-5( VAL ( 2 ) +2*VAL ( 3 > +4*VAL ( 4 ) +2*VAL ( 5 ) + 

X 4*vAL(L)+2*VAL(7)+4*VAL(8>+2*VAL(9)+4*VAL( 10)+2*VAL(1 I )+ 
X4*VAL( 12)+2*VAL( 13)+4*VAL( 14)+2^(-VAL{ 1 5 ) -I-A^VAL ( 16)+2*VAL( 17) 
X +4*VAL( 18)+2*VAL( 19 ) +VAL(20> ) ' ' 

DO 13 J=1 » 10 
PI=J*. I 

9 VAL ( J ) =PHASE ( L ♦ K » J ) -^EXP ( -DELT /P I ) 

13 CONTINUE 

1 1 TEGRAL=« 45/30 0* (VAL( 1 )+4*VAL < 2 ) +2*VAL < 3 ) +4*VAL ( 4 ) 

X +2*VAL(5)+4*VAL(6)+2*VALC7)+4*vAL(S)+2*VAL(9)+VAL( 10) > 

X +TEGRAL 

SOUR { L ♦ K ) =SOUR ( L * K ) +TEGR AL 
RETURN • . - 

END 

©END 
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Chapter VII 


HiOEIDA SYNOPTIC CLIMATOLOGY 


A. OVERVIEW 

The following description of the Climate of Florida has "been 
synthesized from many different published data sources. Special tabula- 
tions of unpublished data have also been constructed to supplement the 
published sets. In particular, extensive use has been made of the tabu- 
lated climatic summaries from Bradley (1972), Newell et al., (1972), 
Baldwin (l97^)» Court (l97^)s various KOAA publications, Air Weather 
Service Climatic Briefs, U.S. Navy Station Climatic Summaries and 
selected NASA Technical Memorandums and Notes. 
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B. SUBFACE SYEpPTIC CLIMATOLOGY 

1. Temperattire 

The normal average daily temperature for January and July is shovn 
in Fig. VII-1 and Fig, VII-2 as taken from Court (l97^)* Note the elimina- 
tion of significant north-south temperature contrast from Janua^ to 
July across Florida. The normal mean daily range of temperature varies 
from 12°C to 10°C north to south across the peninsula in January and 
from 10.5°C to T.5°C in July. 

Figure VII-3 (Court, 19T^) shows the mean interdiurnal variability 
of minimum and maximum temperatures in January and J\ily. In general the 
variability decreases from winter to summer and from north to south as 
continental effects lessen. 

2. Dew Point Temperatoore 

January and July mean dew point temperatures (Court, 19T^) are shown 
in Fig. VII-4 and Fig. VII-5 respectively. Seasonal variations parallel 
those of temperature quite closely. A mean dew point temperature of 
23°C is quite representative of the Florida peninsula in July. 

3. Relative Humidity 

Analysis of relative humidity for 1200 and l800 GMT are presented 
in Fig. VII-6 thru Fig. VII-9 for January and July. These times approxi- 
mate the time of highest and lowest relative humidity respectively. The 
data were gleaned from selected NOAA. and military -climatic summaries with 
variable record periods of at least 10 years. Little difference in the 
1200 GMT values between January and July is noted. A relative inland 
minimum is noted during the early afternoon with slightly drier condi- 
tions prevailing in January. Figure VII-10 taken from Baldwin (19T^) 



VII-3 


shows the hourly relative humidity variation for selected stations based 
upon data for the 195^-1968 period, 

h. Fog 

Days of heavy fog as defined by the occurrence of a visibility 
less than or equal to 0,4 km are shown in Fig, VII-11 and Fig. VII-12 for 
the May-October and November-April periods. The data source is the same 
as lised for the fog frequency maps. Fog is rarely observed along the 
middle and lower east coast and is a maximum in the northwestern part 
of the state in winter. The gradient between coastal and inland loca- 
tions along the peninsula is larger than the figures suggest as fog is 
quite rare along the innaediate coastal strip. 

5. Visibility 

Figure VII-13 end Fig. VII-14 taken from Eldridge (1966) show visibility 
frequencies over the cotenninus -U.S , by season. The data are based 
upon the 1948-1958 period with only Mami, Tampa and Jacksonville data 
used to construct the Florida analyses. The results, however, are 
compatible with synoptic experiences- in that the best -visibilities are 
found along the Florida lower east coast and with the highest values 
found in summer. 

6 . Cloud Cover 

Mean cloud cover from sunrise to sunset for the May-October and 
November-April periods is shown in Fig. VII-15 and Fig. VII-I6. These 
six month periods comprise the approximate wet and dry seasons across 
the Florida peninsula. The data source is the same as previously indi- 
cated with more attention here devoted to small scale variations. 
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A warm season maximum is located along the sea breeze convergence 
region just inland from the lower east coast with a relative minimum 
in the Lake Okeechobee area. Duiring the cool or relatively dry season 
a relative cloud minimum is found in the southwestern part of the penin- 
sula. Daytime cloudiness is a maximum along most of the east coast and 
in the extreme northwest. 

T. Precipitation 

Eainfall maps were constructed from all available civilian and 
military data for the months of March, June, September and December 
with the results shown in Fig. YII-l? thru Fig. VII-20. December 
rainfall is relatively sparse with a gradual decrease from northwest 
to southeast. An exception is the relative maximum along the lower east 
coast. March is the wettest month during the dry season everywhere ex- 
cept along the lower east and southwest coasts and the Keys. A relative 
rainfall maximum is noted from just north of Tampa across the peninsula 
to Cape Canaveral. 

By June the rainy season is well underway with the maximum rainfall 
totals now foimd inland in the lower part of the peninsula. A rather 
large rainfall gradient is noted along the entire east coast (and to 
some extent the west coast) with a relative minimum in the Lake Okeechobee 
region. Convective rainfall is enhanced in the sea breeze convergence 
region approximately 10-20 km inland. The relatively cool waters of Lake 
Okeechobee act to suppress convection, September marks the culmination of 
the rainy season in the southern and western regions of the peninsula 
whereas elsewhere things begin to dry out somewhat. The central and lower 
east coast inland rainfall maximum is still evident as it persists throughout 


the summer 
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The number of days of measureable rainfall (> 0.25 mm) is shown in 

Fig. VII-21 and Fig. VII-22 for the May-October and ITovember-April periods. 

Approximately 50^ of the days record measureable rain during the summer 

along the lower east coast sea breeze convergence zone. This number 

1 

decreases rapidly right along the immediate coast with a value of 36^ 

in the Cape Canaveral area. During the cool, season the inland maximum is 

> 

replaced by a coastal and panhandle maximum. Finally, on an annual basis 
the standard deviation of the number of days with measureable rain in- 
creases from 12 to l6 north to south across the peninsula. 

The mean annual number of days with precipitation > 12.7 ima and the 
mean annual precipitation from days with 12.7 ™i or more as a percent of - 
the total rainfall (Court, 197^) shown in Fig. VII-23 and Fig. VII-2U. 
Comparison of Fig. VII-21 thru lig.. VII-24 suggests that rainfalls reach 
the 12.7 category on 25-30^ of all rain days and account for 70-80^ 
of the annual rainfall total across the peninsula. 

Table VII-I taken from Wallace (1975) shows the diurnal frequency of 
precipitation by amount and category by season. The data soiurce was 
"Climatography of the United States", Wo. 82, "Decennial Census of the 
United States Climate" for the years 1951-1960; published by the U.S. 
Weather Bureau, I962-I963 in addition to Part A of the Uniform Summary of 
Surface Weather Observations, on file at the U.S. Air Force, Environmental 
Technical Applications Center, The data were harmonically analyzed to 
obtain amplitudes and phases of the diurnal and semidiurnal cycles. 
Amplitudes were normalized by dividing them by the 2k hour mean of the 
parameter in q.uestion. 
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Approximately 10-12^ of the summer hours are associated -with pre- 
cipitation yersus J-10% in the ’winter. Trace amoionts occvir h-6% of the 
time whereas heavier amoimts ( >2.5 mm h- occur 2.% and J.% of the time 
in Slammer and ■vri.nter. The time of maximum occtirrence is generally mid- 
afteamoon with the exception of Jacksonville in winter. The noamialized 
an^jlitudes suggest that the afternoon maximum is much more pronounced 
in saammer, Noarmallzed amplitudes in excess of aanity can occaar, e.g, , if 
the precipitation frequency was near zero during the night and high 
daaring the daytime, with symmetary about a noontime peak. 

8 . Thaanderstoami Frequency 

The mean Florida thaanderstoami frequency for the May-Octoher and 
November-March periods is shown in Fig. VII-25 and Fig. VII-26. Civilian 
and aailitary data soaarces were used to constanict the analyses. Additional 
details ad.th special eirrohasis on Cape Canaveral can be found in reports 
by ffeiamann (1968, 1970). The analyses suggest that during the warm season 
thunderstorms are suppressed along the immediate coast and in the Lake 
Okeechobee area, consistent with monthly rainfall statistics. Winter 
frequencies are considerably reduced with the maximum shifting to the 
panhandle area. 

Table ’Vll-rll gives the diurnal cycle in the thunderstorm frequency 
by season for a number of Florida stations based upon data derived from 
Wallace (1975). The percentage of summer hours with audible thunder 
ranges from 3 to 6^ with a maximum across central Florida. A’ very strong 
late afternoon peak is noted for all stations. Winter frequencies are 
considerably reduced across the peninsula (several stations have been 
combined to achieve a reasonable sample size) with a much weaker early 



vii-7 


evening maximinn. A weak morning maximum is noted across the panhandle 
region in winter. Caution must he ^lsed in interpreting this data as 
Wallace notes due to variable record lengths and differing decibel 
levels required for thunder to be heard at the various stations. 

9. Winds 

Prevailing wind directions and mean wind speeds for all available 
civilian and military locations in the state of Florida are shown in Fig. 
VII-27 thru Fig. VII-30 for the months of March, June, September and 
December. Wind speeds average 4 ms~l except near 5 ms“l along portions 
of the lower east coast in March. Directional consistency is weak. By 
J\me the mean wind speeds are reduced to 3-4 ms“l with easterly flow south 
of 28®W and a tendency for southwesterly flow fuirther to the north, 
suggesting an east-west surface ridge line oriented across the mid peninsula. 
By September siorface winds are easterly everywhere with average speeds 
of 4 ms~l. Northerly components giving way to northeasterly components 
from north to south prevail in December with mean speeds in the 4-5 ms“l 
range. Note, the listed wind directions refer to ' prevailing wind 
directions. Lintle information is available on surface wind variability 
for most of the stations listed in the figures. Along the east coast 
wind speeds increase by about 1 ms“l with easterly flow and decrease 
by nearly the same amount with southwesterly flow. For reference surface 
wind roses taken from Baldwin (1974) for January and July are presented 
in Fig. VII-31 and Fig. VII-32. 

10. Fronts 

Figure VII-33 thru Fig. VII-38 show the frequency of cold, warm 
and stationary fronts, in days per 10 years, for the conterminus 
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United States "based upon tabulations from the 1200 GMT Daily Weather 
Map for the 1951-1970 period. The data were presented by Morgan et al. , 
( 1975 ). A north-south gradient of cold frontal passage is quite evident 
with the exception of December. Warm frontal passages across the 
peninsula are con^iaratively rare, Bosart and Korty (1976)^ have provided 
a tentative hypothesis for this observation in terms of the Carolina 
coastal baroclinic zone. Stationary front frequencies suggest a maximum 
across central Florida in all seasons except late spring and su mm er. 

This- is related to the well known tendency for cold fronts to stagnate 
across central Florida. 

11. Stagnating Anticyclones 

Korshover ( 197 ^) has presented a climatology of stagnating anti- 
cyclones for the 1936-1975 period based upon tabxilations from the Daily 
Weather Map published by WOAA. Figures VII-39 thru Fig. VII-i}2 show 
the number of cases of atmospheric stagnation (four days or more) for 
this it-0 year period by 'season. Stagnation is defined for those anti- 
cyclonic regions where surface geostrophic winds are under 7 ms“^ and 
hence actual winds are considerably less. The highest frequencies are 
found in autumn across the southeastern states. Some stagnation cases 
are found across northern Florida in all seasons with a pronoimced 
spring raaximxm across east central Florida. 

12. I Cape Canaveral Climatological Data 

Figure VII-U 3 shows the Air Weather Ser-vice Climatic Brief for 
Cape Canaveral for the 1950-1968 period. This is provided as a 
reference for more detailed information. 

^Published in preprint volume of the Conference on Coastal Meteorology, 
September 1976, Virginia Beach, sponsored by the American Meteorological 
Societv. 
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REPEODUCIBILrrY 

C. UPPEK AIR CLIMA.TOLOGY - FLORIDA ORIGINAL PAGE ^ POOR 

1. Temperature 

Hewell et ai. , (l9T2) have published the only infomiation which 
gives upper air conditions across Florida in which all six radiosonde 
stations are used. The period of record is 1957-196^1 and is based 
mainly on the 0000 GMT observations. 

Table VTI-III gives the 85 O mb, 7 OO mb and 500 mb temperatures by 


winter and simmer season across Florida. The summer uniformity is rather 
remarkable such that detailed Cape Canaveral data should be representa- 
tive of the entire peninsula. As would be expected from synoptic con- 
siderations the coldest winter temperatures are found across the north- 
western panhandle at Eglin AFB (Valparaiso) with a modest thermal 
gradient across the peninsula at all levels. 

Upper air summaries prepared by the National Climatic Center (NOAA) 
for a few of these stations suggest the following standard deviations of 
temperature north to south across the peninsula: (a) 85 O mb: 4.5 to 3.0°C 
in winter and 1.3 to 0.9°C in summer, (b) TOO mb: 3.6 to 2.6°C in 
winter and 1.3 to 0.9‘^C in summer, (c) 500 mb: 3.0 to 2.0*^C in winter 

and about 1.0°C in summer. 

2. Moisture 

Moisture data is not as definitive or as widely published as 
temperature data. Tabulated summaries (0000 GMT data) available from the 
National Climatic Center suggest -the following patterns: (a) 85 O mb: 

60-75^ across most of the Florida peninsula and 40-50^ across the north- 
west panhandle in winter, 70 - 75 ^ and rather uniform in summer, (b) 700 mb: 
25 - 35 ^ from north to south in winter and a fairly imiform 55 - 65 ^ in 
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sxumner, (c) 500 mb: uniform 20-30^ in winter and 4o-50^ in summer. 

Specifics for Cape Canaveral will be presented in a later section. 

3 . Winds 

Table VII-IV gives the u and v wind components (ms“^) at the sirr- 
face, 850 mb, 700 mb and 500 mb for the six Florida radiosonde stations 
based upon the Newell et al. , (1972) data. In addition, standard de- 
viations of the u and v components are given at 850 mb and 500 mb. The 
basic winter west-southwest flow increasing with height is evident above 
the surface with a strong gradient across the Florida peninsula. At 
the surface the picture is more confused but a winter westerly flow across 
the western panhandle giving way to a northerly flow along the' upper east 
coast and a east -northeasterly along the lower east coast is suggested. 

The summer surface picture is probably confused by sea breeze effects 
with an easterly and westerly flow along the east and west coasts 
respectively and a southerly flow along the panhandle area. Aloft an 
east-southeast flow across the lower peninsula gives way to a westerly 
flow across northern Florida. Clearly 3^ east-west oriented ridge line 
, across central Florida in the mean data in summer. This ridge 

line can oscillate north and south as judged by the greater variability 
of the u as opposed to v components in summer at the selected levels. 

Note also the southward slope with height of the boundary between the 
summer easterlies and westerlies. 

4 . Morning and Afternoon Mixing Heights 


Holaworth (1972) compiled data on seasonal and annual mixing 
heights (morning and afternoon) for selected coterminus United States 
stations for the 196O-I964 period. Table 'VII-X shows the summarized 
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data for Tampa, Miami and Jacksonville. No data is readily available 
for the other Florida stations. A general increase in mixing height^ 
is noted from ■winter to summer '^rLth reduced 'wind speeds as convective 
heating becomes more’ dominant . Inland locations would probably have 
lower morning mij^ng heights and somewhat higher afternoon values . 

Table VII-XI summarizes the episode-days of high meteorological 
air pollution potential for episodes lasting at least two days for these 
stations. For example, both Tampa and Jacksonville have 29 episode-days 
in winter of mixing heights imder 1000 m and winds under 6 ms""^ , or 
about 6.5^ of the winter days. The res-ults suggest that the upper 
and lower Florida peninsula is most sensitive to air stagnation in 
winter and summer respectively. 
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D. UPPER AIR CLIMATOLOGY - CAPE GAUATERAL 

1. Temperatures Moisture and Winds 

Table VII-V gires a 20 year summary of temperature, standard de- 
viation of temperature and relative h\imidity for selected lower tropo- 
spheric pressure levels for the indicated months. Likewise, Table VII-VI 
summarizes the data by height levels but only for a 10 year period and 
for different months. As an approximate conversion factor, note that 
the 850 mb and 700 mb levels closely approximate 1,5 km and 3.1 hm 
respectively at the latitude of Florida. Reasonable consistency for 
the Cape Canaveral data is seen amoung Table VII-III, Table YII-V, and 
Table VII-VI despite the varying record lengths. Note especially .the very 
moist conditions in the late summer and early fall. This is typical 
of the east coast of Florida at that time. 

Table VII-VII summarizes Cape Canaveral wind statistics by level 
in more detail based upon a 50^ cumulative frequency of u and V compo- 
nents'. Considerable seasonal variation is evident in the lower tropo- 
sphere. note especially the deep mean easterly flow of September in 
nerms of the potential for the shuttle exhaust cloud to remain over land. 

2, Diurnal Variations 

In an effort to assess the possible significance of diurnal varia- 
tions , ten years of 0000 GMT and 1200 GMT Cape Canaveral data were 
assembled on WBAIT-33 forms provided by the national Climatic Center. 

Table VII-VIII shows the results for the lower troposphere. The 
1200 GMT sounding is uniformly cooler and has a higher relative humidity 
through 900 mb with the differences becoming indistinct by the 850 mb 
level. Little difference in the height field is noted. 



VII-13 


Diurnal wind frequencies by quadrant for June and December are 
shown in Table VII-IX for the 85O mb surface. The results are rather 
interesting and warrant future inwestigation. Differences between on- 
shore and offshore wind components are apparently rather negligible by 
time or by season. However, the 1200 GJ'fT observations feature en- 
hanced southerly flow in comparison to the 0000 GMT observation, 
especially in summer. More months and data for the 0600 GMT and I800 GMT 
radiosonde release times need to be investigated to either support or 
refute this finding. A possible explanation mght be based on differen- 
tial heating whereby the peninsula cools more than the adjacent water 
at night and heats up more during the daytime. This would lead to 
relatively higher and lower 850' mb heights over the inland peninsula 
by late afternoon and daybreak respectively and thus favoring enhanced 
southerly flow along the east coast in the early morning. If this 
explanation has any merit, then Tampa, located along the west coast, 
should have a greater relative northerly flow at the 1200 GMT obsemration 


time. 
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E, SEA. BEEEZE AID RELATED EFFECTS 

During most of the ■warm season the Florida peninsula is relatively 
unaffected hy large-scale organized synoptic scale dist'urbances , 

Weather during this period is then determined hy differential heating 
heWeen land and sea and the attendant interaction with the synoptic 
scale flow. This means that the sea breeze, an almost daily occurrence, 
will play a major role in determining the locai weather patterns. 
Climatological information on the strength, duration, inland and verti- 
cal extent of the sea breeze circulation is lacking. However, the many 
observational and theoretical papers in the literature on the subject 
suggest that the average inland extent of the sea breeze is 10-15 km 
with the vertical circulation confined to the lowest 500-1000 m. 

The picture in Florida is complicated by the existence of a double 
sea-breeze superimposed on small but significant variations in the 
prevailing synoptic flow. Frank et al. , (I 967 ), Frank and Smith (1968) 
and Pielke (1973) have shown that preferred shower locations over 
central and south Florida depend strongly on the location of the sea 
breeze convergence zone. Eecall the inland precipitation ■maximum along 
the entire east coast dtiring the warm season that was discussed pre- 
■viously. On an annual basis Miami International Airport and Miami Beach 
receive l4S and 109 cm of rain respectively, although the distance be- 
tween the stations is only 10 km. The observed difference is almost 
entirely accounted for by warm season precipitation. 

In general on a typical summer day lines of showers form along the 
■windward and leeward coasts of south and central Florida. The windward 
shower line tends to move toward the lee coast during the day and has a 
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tendency to merge ■with the leevard line in the late afternoon, This 
process is particularly evident on days "with moderate easterly synoptic 
scale flow such that extensive precipitation occurs along the "west coast 
sea hreeze convergence. Immediate east coastal locations do not have 
much rain in such, situations. Under light wind conditions, the shower 
lines show little movement during the day, a not too pleasant a prospect 
for those individuals living along the sea breeze convergence zone. 

With a westerly or southwesterly synoptic scale flow the sea breeze con- 
vergence zone is closer to the east coast so that rainfall probabilities 
along the immediate shoreline are enhanced compared to further inland. 
These factors are taken into accoimt in Cape Canaveral thunderstorm 
forecasting and are important to consider when assessing the probability 
of the space shuttle exhaust cloud encountering an active thxmderstorm 
cell. 

Neumann (1968, 1970) has carried out the most definitive study of 
frequency and duration of Cape Canaveral th'understorms . For reference 
pirrposes Fig. VTI-Ui^ shows a map of the Cape Canaveral area. A similar 
s'ummer precipitation gradient prevails between Titusville and the ocean 
front, although somewhat weaker, as compared to Miami Airport and 
Miami Beach. 

Figure 711-45 and Fig. 711-46 show the prominent daily summer 
thunderstorm peak reaching 50^ in July and August in the vicinity of 
the Kennedy Space Center. The strongly diurnal nature of the percent 
of days with a thunderstorm is quite evident. The data in these figures 
is based upon 13 years: 1951-1952 and 1957-1967. It is also instruc- 

tive that Neumann finds that 11^ of the thunderstorm occurences at 
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Cape Canaveral take place entirely between regular hourly surface re- 
ports. This implies a significant underestimation of thunderstorm fre- 
quencies in the published climatological summaries (e.g. , ¥ad.lace [1975]) 
which are based entirely on hoiorly observations. 

Figure YII-4 y shows that the beginning and end of the Cape Canaveral 
thunderstorm season is associated with a strengthening and weakening of 
the 1200 GMT 1000 m resultant wind, consistent with the ideas expressed 
previously. More details are provided in Fig. VII-48 on the annual 
variation in the 1200 GMT 1000 m vector and scalar winds imder conditions 
with and without afternoon thunderstorms . Afternoon thimderstorms are 
clearly favored by westerly wind components and suppressed by easterly 
components during the summer. Afternoon thunderstorms during the cool 
season are especially prevalent in the strong southwesterly flow ahead 
of the more potent migrating troughs at that time of the year. 

The data is summarized in Pig, VII-49 which shows the probability 
of afternoon thunderstorms over the entire May-September thunderstorm 
season as a function of the 1200 GMT 1000 m wind speed and direction. 

The contrast between northeasterly and southwesterly flow is remarkable. 
Southwesterly flow days are associated with a sea breeze convergence 
zone very close to the coast and the developing thunderstorm cells then 
can drift towards the ocean during the afternoon. Finally, Fig. VII-50 
shows a forecast monogram for the peak of the thunderstorm season 
around 1 August. ¥ith a wind of 250 degrees between 5 and 9 ms“^ the 
afternoon thunderstorm probability reaches 
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F. COHCLUSIOHS 

Tlie data presented in the previous sections can be used for many 
individual purposes. In terms of the space shuttle program it becomes 
impornant to worry about possible inadvertent weather modification by 
the stabilized ground cloud. Clearly, the risk is greatest when the 
exhaust cloud remains over or near land, especially populated areas. 
These risks would be enhanced under the following synoptic conditions: 

1. Surface sea breeze regime and weak westerly flow at 1000 
and 2000 m in summer 

2. Pre cold front squall lines that occasionally occur in the 
winter and early spring 

3. Post cold frontal (weak with easterly surface flow giving 
way to weak westerly flow aloft 

U. Pre warm front passage in winter with low land onshore flow 
veering to offshore flow aloft 

5. Weak stationary front located south of Cape Kennedy with sur- 
face onshore flow veering to offshore flow aloft 

6. Stagnating anticyclonic conditions 

7. General easterly flow of late summer and early fall across 
the Florida peninsula 

8. Thunderstorm events 

9. Tropical storm situations 

Items (l) through (3) above are in particular quite difficult to 
evaluate without further research. The- surface wind field on sea 
breeze days is relatively well known. Unfortunately, the lower tropos- 
pheric wind field is not at all well documented in terms of diurnal and 
seasonal variations as a function of the prevailing synoptic patterns. 
Consequently, accurate mesoscale three-dimentional trajectories are 
rather difficult to compute. 
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Pre cold front sguall lines may occasionally move across northern 
and central Florida in the 'vTinter and early spring. These sguall lines 
are usually oriented from northeast to southwest and may propagate 
southeastward at speeds up to knots . The freguency of occurrence 
may reach two to four times per month at the height of the cool season. 
Such sguall lines usually, die out before they reach south Florida. 

Strong low level southerly flow ahead of the sguall lines carries the 
potential risk of incorporating particulate matter from the space 
shuttle exhaust into the active cumulonimbus clouds of these systems. 

The risk would be greatest within the- first 2i(— hours after launch. 

The results from Morgan et al. , (1975) suggest possibly two warm 
front egisode days per winter month across central and northern Florida. 
Assuming a 12-hovir frontal passage this would work out to be about a maxi- 
mum warm front influence of 2>% during a mid-winter month. Convective 
activity often accompanies warm frontal passage especially in northern 
Florida. Pre frontal veering winds can carry the exhaust cloud north- 
ward into the region of favorable convection. 

Stationary fronts are more ambiguous because the lower tropospheric 
wind field can be considerably different depending upon such factors 
as frontal intensity and distance from the station. An upper bound of 
4 days per winter month gives an upper bound of a 13^ influence and 
this will be an overestimate of the true risk. 

The results from Korshover (1976) and Holzworth suggest that re- 
latively stagnant conditions are possible in winter and early spring 
across central and northern Florida from 3 to of the time, ranging 
down to 1 to 2.% in summer, formally such situations are associated with 
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stable air but occasional air mass thunderstorms during the warm season 
pose' some risk for vertical mixing. 

The deep easterly flow regime of late summer and early faJ.1 becomes 
a special problem if the exhaust cloud encounters active precipitation 
cells after launch. This risk is difficult to assess with the data 
currently in hand. Overall, the thunderstorm problem is perhaps 
potentially most serious because of uncertainties in precise prediction. 
Neumann (1968, 19T0) has compiled a very valuable data set which 
obviously is heavily used in thunderstorm prediction at Cape Canaveral. 
His data suggest that 1200 CSIT 1000 m level winds of 1 ms“^ or less 
can be encoiuitered of the time during the May-September period. 

This number increases to 15.1^ for a 2 ms“^ or less cut off point. This 
suggests in t^^rn that the exhaust cloud can still be within 10 km of 
Cape Canaveral more than an hour after launch, assuming relatively 
light surface winds, on a substantial number of simmer mornings. How 
this simple calculation would change based upon the 0000 GMT data is 
unknown now and is worthy of further research. 

Finally, tropical storm situations carry potential risk because 
of the active cumulonimbus clouds associated with such systems and 

I 

potential for inadvertent weather data. Historical data (Court, 197^) 
suggest that from 10-15^ of the June-October Florida rainfall is 
associated with tropical disturbances. 
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G. KECOMMENDATIOEfS FOR FURTHER RESEARCH 

The most pressing need is to establish typical and atypical three- 
dimensional air trajectories in the lower troposphere as’ a function of 
the time of the year and prevailing synoptic situation over the Florida 
peninsula. Likewise, temporal and spatial variations of various 
meteorological parameters on scales of several hours and several tens 
of kilometers needs to he better documented. In particular, the following 
is suggested. 

1. Diurnal and seasonal variations of the lower tropospheric wind 
field across Florida using the 0000, 0600, 1200 and 1800 GMT 
winds aloft observations. Published studies rarely refer 

to winds other than at the 0000 GMT observation time.' 

2. Diurnal and seasonal variations of various meteorological 
surface parameters along a line perpendicular to the east 
coast to better understand the complex mesoscale air motions. 
Existing stations in the Jacksonville area offer such an 
opportunity for research. The present upper air network does 
not permit such a study above the surface layer, however. 

3. Establish hail frequencies across Florida on a county by 
county basis to better assess the risk of inadvertent weather 
modification diiring convective situations. 

Better establish precipitation frequency, intensity and 
duration as a function of the time of day and season from the 
existing hourly precipitation network across Florida. 
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Figure VII-2; July Mean Temperature 
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Fig.37. Mean interdiurnal variability of minimum and maximum temperatures ('Q in January and July, based on data at 75 stations, 1957-1961. (Prepared by 
Stanley Marsh from data of Landsberg, 1966 .) 


Figtire VII-3: Mean Interdiurnal Variability of Minimum and Maximum Temperatures (°C) 
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Climate of the conterminous United States 





Fig.lO, Mean monthly dew point (*C) irt January (A) and July (B), based on observations every 
V. Dodd’ from (!ata used for Dodd, 1965.) 


Figure VII-5: 
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5’igui‘e VII— 6; Surface Mean % Relative Humid.ity 
January 1200 GMT 












Figure VII-7: 


Surface Mean % Relative Humidity 
January 1800 GMT 
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Figure VII-8: 


Surface Mean % Relative Huraidity 
July 1200 GMT 










Figure VII-9: Surface Mean % Eelative Humidity 

July 1800 GMT 
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Figure VII-10: Mean Hourly % Relative Humidity; Jan. (dashed) 

7 July (solid) 
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Figure VII-11: 


Days of Heavy Fog: 
May-OctoHer 


Visibility < O.U km 















Fig. 3. Percentage of time by seasons the visibility is less than 5.0 km. » 

Figure VII- 13 : Surface Visibility- 

After Eldridge (1966) 












FlC. 4. Percentage of time liy seasons tlic visi!)ilily is less than 10 knj. 



Fig. S. Percentage of time by seasons the visibility is less than 20 kia. 


Figure VII-lU : Surface Visibility 


After Eldridge (1966) 











Figure VII-15: Mean Cloud Cover (Tenths) Sunrise -Sunset 

May-Octoher 
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Figure VII-17: Mean Monthly Precipitation (mm) 

December 
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Figure VII-I8: 


Mean Monthly Precipitation (mm) 
March 










Figure VII-19: Mean Monthly Precipitation (sun) 

June 
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Figure VII-21: Mean Number of Days Precipitation > 0.25 nmi 

May-October 








Figure VII-22: Mean Number of Days Precipitation > 0.25 

N ovember-Apri 1 














Figure VII-23: Mean Annual Number of Days with Precipitation of or more 
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Fig.l6. Mean annual number of days with precipitation of 12.7 mm or more, based on data for 315 stations, 1942-1961 (A), and ratio (%) of total annual rainfall i 
such days to total from all days (B) based on data for3I5 stations, 1942-1961. (From I'avlhus and Miller, 1964.) 


Figure VII-2^r Mean Annual Precipitation from Pays 'with. 12.7 nun or more as percent of Total Amount 
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Figure VII 
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Figure VII-28; 
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Figure VI 1-29: 


Prevailing Surface Wind Direction and Speed (ms 
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Figure VII-30: 


Prevailing Surface Wind Direction and Speed (ms-^) 
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Figure VII-31: Mean Surface Wind Roses: January 


After Baldwin (l97^) 













Figure VII-32; Mean Surface Wind Roses: July 
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Figure VII-33: Cold Front Climatology 

■ March j June 


After Morgan et al. (1975) 





Figure VII-3^: Cold Front Climatology 

September & December 


After Morgan et al. (1975) 




Figrire VII-355 Warm Front Climatology 
March & June 


After Morgan et al. (1975) 



Figure VII-36: Warm Front Climatology 

September & December 
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Figure VII-37! Stationary Front Climatology 


After Morgan et al. (19T5) 




Figirre VII-38; Stationary Front Climatology 
September & December 


After Morgan et al. (1975) 
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Figure VII-UO: Sprint Atmospheric Stagnation 
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Figure VII-4l: Summer Atmospheric Stagnation 
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Figure VII-U2: Autumn Atmospheric Stagnation 
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Figure VII~43; Cape Canaveral Climatic Summary 
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Figure L3. Map of Cape Kennedy a,rea. 


Figure VII-44: Cape Canaveral Area Map 
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Figure la. Probability (%) of at 
least one thunderstorm on August 1 g 
(EST) between time and time Tq | 

+ A T. (data derived from Part I). I 

£ 
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Figure lb. Daily thunderstorm 
frequencies (top panel) smoothed 
over periods of 5, 15, and 31 days 
(data derived from Part 1). 

After Neumann (1970) 
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Figure VII-45: Smoothed Cape Canaveral Daily Thunderstorm Frequencies 

(lover Panel) 
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Figure 2. Probability of thvinderstorms at or in the immediate vicinity of 
'the Kennedy Space Center over specified time intervals (data derived from 
Part I). 


Figure VII-46: Thunderstorm Probability at Cape Canaveral as a Function of Time 












































After Uemaann (1970 ) 


Figure 4. Location ( ^) of the 1200 GMT 3000-foot resultant wind at the 
Kennedy Space Center for each of the 73 dates referred to in figure 3. The 
location ( © ) of the resultant wind for the 15th day of each month is inter- 
polated from the location of the adjacent 5-day positions. The location ) 
of the resultant wind for the entire thunderstorm season is 187 degrees at 
3. 6 knots. 


Figure VII-47: Cape Canaveral 1200 G-MT 1000 meter Resultant Winds during the Thunderstorm Season 
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After Neumann (l970 ) 

Figure 6, Annual variation in the 1200 GMT 3000 -foot vector and scalar 
winds under conditions with and conditions with and without afternoon thunder- 
storms. The vector locations of the wind components under conditions with- 
out afternoon thunderstorms are given for the period' May through September 

only. 

Figxn'e VII-48: Cape Canaveral 1200 GMT 1000-m Vector and Scalar Wind Variations 

During the Summer Season. 
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Figure VII-50: Sampel Cape Canaveral Thunderstorm Forecasting Monogram 
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TABLE VII-I (Wallace - 1975) 


Diiornal Frequency of Precipitation by Amount and Category by Season 



Percentage of hoirrs 






with Precipitation 

Phase of Diurnal Cycle 



Station 

(Trace or More) 

(Time of Maximum — GMT) 

Normalized 

Amplitude 


Jun-Aug 

Hov-Mar 

Jun-Aug 

Hov-Mar 

Jiin-Aug 

Nov-Mar 

Jacksonville 

10 

10 

2100 

1300 

.90 

.10 

Orlando 

12 

9 

2100 

2000 

1.10 

.15 

West Palm Beach 

11 

9 

1800 

2000 

.65 

.20 

Miami 

10 

7 

1900 

1800 

.40 

.10 


Percentage of Hours with 



It 



Precip > 

2.5 mmh“'‘ 


It 

IF 



■Jun-Aug 

Kov-Mar 

Jun-Aug 

Wov-Mar 

Jun-Aug 

Nov-Mar 

Jacksonville 

1.8 

1.1 

2100 

l4oo 

1.05 

.20 

Orlando 

2.2 

1.1 

2100 

2300 

1.15 

.40 

West Palm Beach 

2.1 

1.1 

1900 

2100 ■ 

.70 

.35 

Miami 

2.0 

0.9 

2100 

2000 

.60 

.25 

■ 

Percentage 

of Hours with 






a Trace 


IT 

II 



Jun-Aug 

Nov-Mar 

i 

Jun-Aug 

Kfov-Mar 

Jun-Aug 

'Nov-Mar 

Jacksonville 

5 

5 

2100 

1300 

.60 

.05 

Orlando 

6 

4 

2100 

2100 

.55 

■ .20 

West Palm Beach 

5 

5 

1800 

2000 

.50 

.15 

Miami 

5 

4 

1900 

1800 

.35 

.20 



TABLE VII-II (Wallace - 1975) 

Diijmal Cycle in the Thunderstom Frequency by Season 



Percentage 

of Hours 

Station 

Audible 

Thunder 


Jun-Aug 

Wov-Mar 

Eglin AFB 

U.5 

5.5 

Panama City 

3.P. 

4.8 

Tallahassee 

k.8 

7.6 

Jacksonville 

2.8 

1.7 

Daytona Beach 

4.9 ') 

1 

3.2 

Orlando 

3.8 J 

1 

Tampa 

4.4 

1 

) 

Vero Beach 

5.6 ] 

f 2.2 

Fort WJyers 

4.4 X 

) 

West Palm Beach 

3.2 -n 

. 2.0 

Miami 

3.7 i 

1 


Phase of 

Diurnal 

Cycle (Time of Max-GMT) 

Jun-Aug 

Wov-Mar 

2000 

0300 

2000 

1000 

2100 

1600 

2100 

2300 

2100 1 

t 2300 

2200 / 

1 

2200 "" 

) 

2100 j 

r 

0000 

2200 J 

1 

2100 “1 

} 

► 0000 

2100 J 

[ 


Hormalized Amplitude 
Jim-Aug Kov-Mar 
.65 
.70 
1.30 
l.UO 

1.50 
1.60 
1.55 
1.40 
1.60 
1.10 
1.10 




TABLE VII-III 


Station 

Eglin AFB 

Jacksonville 

Cape Canaveral 

Tampa 

Miami 

Key West 


Mandatory Level Temperatures Across Florida 
After Newell, et al. (l972) 

850 mb 700 nib 500 mb 


Dec-Feb 

Jun-Aug 

Dec-Feb 

Jun-Aug 

Dec-Feb 

Jun-Aug 

7.6 

17.5 

1,5 

8.1 

-14.2- 

-7.1 

8.0 

17.5 

1.8 

8.4 

- 13.9 

- 7.2 

10.0 

17.7 

3.8 

8.1 

-12.3 

-7.lt 

9.6 , 

17.8 

3.6 

7.9 

- 12.0 

-7.3 

11.3 

17.lt 

5.1 

8.2 

-10.4 

-7.lt 

12.2 

17.8 

5.7 

8.4 

1 

CD 

• 

- 7.1 


Units {°C) 



TABLE VII-IV 


Mandatory Level Winds Across Florida After Newell et al. (l9T2) 

Dec erab er-February 


Station Surface 85 O mb 700 mb 500 mb 



u 

V 

u 

o{n) 

V 

a(v) 

u 

V 

u 

a(u) 

V 

a(v) 

Eglin AFB 

0.2 

-0.0 

5.2 

6.9 

0.4 

7.8 

13.0 

1.5 

23.3 

10.8 

3.5 

10.9 

Jacksonville 

0.1 

-0.7 

6.3 

6.6 

1.4 

6.9 

12.8 

1.5 

22.9 

10.9 

2.1 

9.7 

Cape Canaveral 

0.5 

-0.8 

8.1 

6.1 

2*2 

8,7 

10.1 

1.7 

29.4 

6.1 

5.1 

i4.5 

Tampa 

1.0 

-1.2 

3.7 

6.7 

1*1 

5.8 

15.7 

2.8 

18.1 

8.2 

3.6 

7.8 

Miami 

-0.8 

-0.2 

2.2 

6.9 

0.6 

4.9 

7.2 

1.6 

15.5 

8.6 

5.0 

7.0 

Key West 

-1.1 

-1.7 

1.5 

6.5 

0.3 

5.1 

6.1 

1.0 

17.0 

8.3 

1.7 

6.2 


June-August 



u 

V 


a-(u) 

V 

a(v) 

u 

V 

u 

o(u) 

'v 

a(v) 

Eglin AFB 

0.9 

1.8 

0.4 

4.6 

0.2 

4.1 

1.3 

0.0 

2.7 

5.9 

0.1 

5.0 

Jacksonville 

-1.9 

1.4 

1.2 

4.6 

0.5 

4.1 

1.6 

0.3 

2.5 

5.6 

0.2 

4.7 

Cape Canaveral 

-1.8 

0.4 

-1.0 

6.6 

1.7 

5.0 

1.3 

1.4 

2.0 

8.4 

0. 6 

6.0 

Tampa 

0.5 

- 0.3 

- 0.7 

4.4 

1.1 

3.5 

1.0 

0.8 

1.0 

5.4 

0.6 

4.3 

Miami 

-1.'9 

0.9 

-1.6 

4.3 

0.6 

3.2 

- 0.7 

0.9 

- 0.4 

5.2 

1.1 

4.0 

Key West 

-2.0 

0.6 

-2.3 

4.1 

0.7 

3.1 

- 1.0 

0.9 

- 0.9 

4.8 

0.5 

3.5 


UNITS; ms”^ 



TABLE VII-V 


Cape Canaveral; Selected Upper Air Data by Pressure 




Jan 

Apr 

Jul 

Oct 

1000 mb 

T(°C) 

15.7 

20.5 

25.6 

23.2 


(T) 

4.8 

2.9 

1.7 

2.7 


^EH 

73 

70 

80 

75 

850 mb 

T(°C) 

9.3 

12.5 

17.3 

14.3 


(T) 

3.7 

3.1 

1.2 

2.5 


%m 

56 

56 

71 

69 

700 mb 

T(°C) 

3.1 

5.0 

7.9 

6.7 


(T) 

3.1 

2.4 

1.2 

2.1 


^RH 

33 

34 

59 

46 


Data Source: Patrick AFB 1950-1956 

Cape Canaveral 1956-1970 


0000 GMT Only 



reproducibility of the 


TABLE VII-VI 


Cape Canaveral: Selected Upper Air Data by Height 

March June September December 



T(°C) 

^EH 

T(°C) 

^RH 

T(°C) 

^RH 

T(°C) 

^RH 

Surface 

20 . T 

81 

27.5 

79 

27.9 

80 

18.5 

85 

500 m 

'IT. 5 

?4 

23.2 

80 

23.7 

83 

16.1 

78 

1000 m 

i4.7 

75 

20.0 

72 

19.7 

79 

13.2 

78 

1500 m 

12.4 

75 

17.2 

73 

17.2 

78 

■ 10.8 

75 

2000 m 

10,1 

63 

13.9 

69 

i4.7 

74 

9.4 

57 

3000 m 

5.0 

4l 

8.5 

65 

9.4 

70 

5.1 

25 

4000 m 

- 0.6 

33 

2.9 

42 

3.9 

64 

- 0.3 

— 

5000 m 

- 6.9 

33 

- 2.8 

42 

- 1.8 

57 

- 6.3 



Data Source: MSA - Cape Canaveral ( 195 O-I 96 O) 

0000 GMT Data Only 



TABLE VII-VII 


Cape Canaveral: 50^ Ciamulative Frequency Winds by Height 


March 



u 

V 

9 

Vr 

Siirface 

- 0.4 

-0.4 

045 

0. 6 

1000 m 

2.6 

0.8 

253 

2.7 

2000 m 

6.0 

0.6 

264 

6. 0 

3000 m 

9.4 

0.4 

268 

9.4 

UOOO m 

13.5 

-0.3 

271 

13.5 

5000 m 

17.2 

-0.3 

271 

17.2 


September 



u 

V 

0 

V 

''r 

Surface 

- 2.8 

-0.6 

078 

2.9 

1000 m 

- 3.7' 

-0.4 

084 

3.7 

2000 m 

- 2.6 

-0.4 

08l 

2.6 

3000 m 

- 1.3 

-0.3 

077 

1.3 

4000 m 

- 1.8 

-0.2 

084 

1.8 

5000 m 

- 0.7 

-0.4 

060 

0.8 


UNITS; ms“^ 

0 = resultant wind direction 

= resultant wind speed 
|V| = scalar mean wind speed 


June 


|v| 

u 

V 

0 


|V| 

■3.7 

- 0.9 

-0.0 

090 

0.9 

2.8 

7.4 

- 0.3 

0.7 

156 

0.8 

4.6 

8.4 

0.5 

0.0 

270 

0.5 

4.6 

11.7 

1.0 

-0.2 

281 

1.0 

4.5 

15. -5 

1.2 

-0.4 

288 

1.3 

4.6 

19.2 

1.3 

-0.4 

287 

1.4 

4.8 


December 


Iv| 

u 

V 

0 

Vr 

|vi 

3.0 

- 0.0 

-1.4 

360 

1.4 

3.1 

5.5 

- 0.6 

-0.7 

o4i 

0.9 

7.1 

5.2 

1.7 

-0.2 

083 

1.7 

6.3 

4.8 

5.0 

-0.0 

270 

5.0 

8.1 

4.8 

9.1 

0.2 

269 

9.1 

11.4 

4.9 

11.8 

0.4 

268 

11.8 

i4.3 


Data Soiorce; 1950-1960 - NASA 
0000 GMT Only 



TABLE VII-VIII 


Cape Canaveral: Diurnal Variation of Height 

Temperature and Relative Humidity- 


Time 


GMT 

Surface 
PEES (mb) T(°C) 

^RH 

HT(m) 

950 mb 

t(°c) 

^RH 

HT(m) 

900 mb 
T(°C) 

^RH 

HT(m) 

850 mb 
t(°c) 

^RH 

March 

0000 

1017.0 

18.4 

73 

587 

16.1 

65 

1045 

13.1 

66 

1524 

10.3 

63 

1200 

1017. 7 

l4.8 

85 

588 

l4.6 

71 

1045 

12.1 

67 

1522 

10.0 

60 

J-une 

0000 

10ll;.6 

25.9 

78 

583 

22.3 

74 

1052 

19.6 

71 

1542 

16.6 

70 

1200 

1015.2 

23.3 

89 

585 

21.4 

79 

1058 

18.9 

72 

1542 


70 

September 

0000 

1013.7 

26.1 

80 

577 

22.5 

79 

1046 

19.7 

77 

1536 

16.8 

76 

1200 

1014.0 

23.7 

89 

577 

22.0 

83 

io46 

19.2 

79 

1536 

16.5 

76 

December 

0000 

1018.6 

15.9 

79 

596 

14.5 

67 

1052 

12.0 

64 

1529 

9. '7 

59 

1200 

1019.1 

12.9 

85 

595 

13.6 

72 

1049 

11.6 

66 

1528 

9.7 

58 


DATA SOURCE: WEAK 33 Poms I96O-I969 

National Climatic Center 



TABLE VII-IX 


Cape Canaveral; 85O mb Diurnal Wind Frequencies (%) 


Time (GMT) 


Direction Frequency {%) 


0 ) 

,8 

CO 

d 

o 


(1) 

o 

,d 

ca 

O 


Xi 


I 

m 


June 


Dec. 



001-090 

091-180 

181-270 

271-360 


-(%) — 


0000 

23 

2 k 

29 

24 

47 

53 53 

47 

1200 

19 

26 

44 

11 

45 

55 30 

70 


0000 

14 

15 

35 

36 

29 

71 

50 

50 

1200 

12 

17 

4 o 

31 

29 

71 

43 

57 


DATE SOURCE: WEAR 33 Forms 196O-I969 

Rational Climatic Center 



TABLE VII-X 


Florida Mean Seasonal and AnnuaJ. Morning and Afternoon 
Mixing Heights (H) and Wind Speeds (u) for HOP^ and ALL^ Cases After Holzworth (1972) 


! 

Station 

o 

o 

E 

H 

■m 

Winwi 


Spring 

Summer 

Autumn I 

• Annual 


■ 

U, msec ^ 

H, m 


U, m sec“^ 


54 NOP 

U, m scc'*'^ 

m 

1 

U, m sec’“^ 

H.m 

m 

— 1 

U. m soc 

NO? 

All 

Hi 

iiii 

NOP 

All 

NOP AM 

NOP 

All 

NOP 

All 1 


NOP 

AM 

NOP All 

1 

: NOP 

All 1 

NOP 

All 

%N0P 

NOP 

All 




* 























JacksonviUo. 

M 

345 

403 

79.4 

5.2 

5.9 

447 477 

90.4 

5.3 

5.6 

667 

583 

91.1 


4.3 

44 

418 468 

85.9 

4.7 

50 

444 

480 

86.7 

4.9 

5.2 

Florida 

A , 

loss 

1104 

80.1 

8.7 

7.0 

1G39 1667 

86.1 

7.1 

7.2 

1681 

1712 

68.0 


5.6 

58 

1321 1342 

80.4 

6.5 

65 

1424 

1456 

78.6 

6.5 

6.7 


Tampa, 

M 

„394 436 

85.8._ 

-5,8 

—6,1 

503- 526 

91.7“ 

— 5.6- 

— 5.8" 

—656—674 

9Vt“ 


— ^419^ 

4(19 

M 9 





^ Florida 

A 

1052 1079 

81.4 

64 

6.6 

1523 1544 

87.8 

6.7 

6.8 

1460 1526 

63.9 

5.0 5.3 

1401 

1429 

04.4 

6.4 68 

1359 

1394 

80.6 


Miami, 

* 1 
M 

654 

707 

87.2 

5.4 

5.7 

947 

980 

91.1 

5.7 

5.9 

1041 

1071 88.3 

4.3 

4.5 

872 

933 

82.4 

s.o 

53 

878 

923 

! 

87.2 

5.1 

5.3 

Florida 

A 

1208 

1221 

89.2 

0.4 

6.5 

1440 

1459 

S7.4 

68 

6.9 

1360 

1383 73 7 

5.3 

5,5 

1315 

1341 

78.7 

6.6 

69 

1330 

1351 

62.2 

6.3 

65 


DATA SOURCE: 1960-1964 

%0P: Omitting Cold Advection and Precipitation Cases 

^ALL: Includes all Available Data 




























TABLE VII-XI 


Episode-Days of High. Meteorological Air Pollution Potential 

(After Holzvorth [1972]) 



Wind Speed 


Mixing Heights (m) 


Station 

(ms-i) 

500 

10000 

1500 

■ 2000 


2 

0(0) 

0 (0) 

0 (0) 

0 (0) 

Jacksonville 

4 

0(0) 

2 (5) Win. 

12 (28) Win. 

23 (50) Win. 


6 

1(3) Win. 

12(29) Win. 

43(116) Win. 

105(314) Aut. 


2 

0(0) 

0 (0) 

0 (0) 

0 (0) 

Tampa 

4 

0(0) 

-1 (2) Spr. 

9 (23) Win. 

15 (38) Sum. 


6 

1(2) Win. 

14(29) Win. 

79(234) Win. 

127(406) Sum. 


2 

0(0) 

0 (0) 

0 (0) 

0 (0) 

Miami 

4 

0(0) 

0 (0) • 

17 (38) Sum-Aut. 

25 (60) Sum. 


6 

0(0) 

3(11) Win. 

79(207) Sum. 

128(4o 6) Sum. 


DATE SOUECE: 1960-1964 


First figure is the number of episodes j the number of episode-days is given in parentheses. 
Seasonal Peaks as Indicated. 



Chapter VIII 
SUJTHSSIS 


The following risk situations for inadvertent weather modification 
due to the space shuttle exhaust are now sunmarized 

1. Exhaust cloud encountering active convective precipitation cells 
with conseauent vertical transport to the upper troposphere and potential 
for acid rain 

(a) sea breeze convergence during the warm season with attendant 
affcernoon thunderstoxmis . Effects include possible localized 
hail and brief wind gusts in excess of 20 ms-^. Affected area 
is less than 100 km^ with a time scale of less than T + 1 day. 

(b) frontal and prefrontal activity including squall lines mth 
attendant th\mderstorms . Effects include possible localized 
hail, wind gusts in excess of 30 ms“^ and tornadoes. Affected 
area is 100-500 km^ with a time scale of less than T + 2 days. 

(c) general air mass thunderstorms not associated with (a) and (b) 
above but responding to different summer synoptic flow patterns. 
Effects include possible localized hail and brief wind gusts 

in excess of 20 ms~^. Affected area is less than 100 km^ with 
a time scale of less than T + 1 day. 

(d) tropical storms in the vicinity of the Florida peninsula 
within 2k hoxirs of launch- time. Potential effect of shuttle 
exhaust cloud cau^t up in the circulation of a tropical storm 
is unknown in terms of inadvertent weather modification. A sub- 
sequent change of direction of such a storm might be interpreted 
as not an ”act of God" by some people with possible social and 
legal problem from communities in the landfall region. 

lOTE: Approximately 10-12^ of the summer hours are associated with 

precipitation versus 7-10^ in winter across Florida. Trace 
amomits occur k~6% of the hours whereas heavier amounts 
(> 2.5 mmh'’^) occur 2% and 1^ of the time respectively in Florida. 
The time of maximum occurrence is generally midafternoon with 
the exception of northern Florida in winter. The percentage 
of summer hours with audible thunder ranges from 3 to 6^ with 
a maximum across central Florida. A very strong late afternoon 
peak is noted for all stations. Winter frequencies are considerably 
reduced with a much weaker early evening maximum. Overall, these 
numbers should provide an overall upper bound of the percentage 
chance of encountering an active ctimulonimbus cloud at any 
random time. 

2. In the months November-April, when advective and radiative fogs 
maximize, very significant worsening of visibility conditions in foggy 
situations could occur within the area affected by the dissipating S.G.C. up 
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to T + 1 day (area affected up to 10^ 3sm^) and particularly under wind 
flow conditions from the S-S quadrant, 

3, Minor risk associated with easterly flow in lower troposphere 
(unless tropical disturhances are present) particularly in those situations 
where atmosphere is stable under those conditions, clouds do not reach 
the level where ice phase processes are operative. However, overseeding 
of warm clouds with CCK could result in a very significant reduction of 
precipitation over the entire area affected by the dispersing cloud. 

Effect diminishes after T + 1 day. (Criteria; shallow warm cloud 
system and no ice phase. ) 

Stagnating anticyclonic conditions with reduced dispersion of 

S.G.C. Little cloudiness is nonaally associated with conditions of this 
type. The impact is therefore restricted only in the area of visibility 
detoriation and solar energy reduction. This constitutes therefore a 
nuisance and conceivably might violate EPA standards . On rare occasions 
air mass thunderstorms may develop, particularly along the sea breeze 
convergence zone, under stagnant anticyclonic conditions during the 
warm season. The risk would then be equivalent to l(c) above. 

5. Possible modification of a major hurricane located east of 
Florida peninsula at time of launch. Air from latmch site would participate 
in the storm circulation and might indeed cause some modification effects 
producing unknown results. Any subsequent veering of such a storm would 
undoubtedly cause serious social and legal problems. 

6. Cumulative effects; for the projected UO launches per year 
ass Timing several days spacing between launches is considered negligible. 

7. Minimal risk and impact: strong westerly winds system extending 

through the lower troposphere. 
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REASSESSMENT CRITERIA 

A reassessment is required if: 

(a) the volume of the stabilized ground cloud differs by- 
more than 50^ of the value assumed here 

(b) the number of giant particles in the stabilized ground 
cloud with radius larger than 20 y exceeds 1 particle 
per 1000 cm^ 

(c) "the number density of ice nuclei in the stabilized 

ground cloud exceeds 5^ 'of the total number of particles 
present in the S.G.C. (irrespective .of particle size) 

(d) the launch area is moved out of the Florida area 
(our assessments relies heavily upon the climatology 
of the launch area) . 



